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Abstract 

A search for pair production of vector-like quarks, both up-type (T) and down-type (B), as 
well as for four-top-quark production, is presented. The search is based on pp collisions 
at yfs = 8 TeV recorded in 2012 with the ATLAS detector at the CERN Large Hadron 
Collider and corresponding to an integrated luminosity of 20.3 fb” 1 . Data are analysed in the 
lepton-plus-jets final state, characterised by an isolated electron or muon with high transverse 
momentum, large missing transverse momentum and multiple jets. Dedicated analyses are 
performed targeting three cases: a T quark with significant branching ratio to a IT boson and 
a Zr-quark (TT —> Wb+X), and both a T quark and a B quark with significant branching ratio 
to a Higgs boson and a third-generation quark {TT —> Ht+X and BB —> Hb+X respectively). 
No significant excess of events above the Standard Model expectation is observed, and 95% 
CL lower limits are derived on the masses of the vector-like T and B quarks under several 
branching ratio hypotheses assuming contributions from T —> Wb, Zt , Ht and B —> Wt, Zb, 
Hb decays. The 95% CL observed lower limits on the T quark mass range between 715 GeV 
and 950 GeV for all possible values of the branching ratios into the three decay modes, and 
are the most stringent constraints to date. Additionally, the most restrictive upper bounds on 
four-top-quark production are set in a number of new physics scenarios. 
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1. Introduction 

The discovery of a new particle consistent with the Standard Model (SM) Higgs boson by the ATLAS [1] 
and CMS [2] collaborations is a major milestone in high-energy physics. However, the underlying nature 
of electroweak symmetry breaking remains unknown. Naturalness arguments [3] require that quadratic 
divergences that arise from radiative corrections to the Higgs boson mass must be cancelled by some 
new mechanism in order to avoid fine-tuning. To that effect, several explanations have been proposed in 
theories beyond the SM (BSM). In supersymmetry, the cancellation comes from assigning superpartners 
to the SM bosons and fermions. Alternatively, Little Higgs [4, 5] and Composite Higgs [6, 7] models 
introduce a spontaneously broken global symmetry, with the Higgs boson emerging as a pseudo-Nambu- 
Goldstone boson [8]. Such models predict the existence of vector-like quarks, defined as colour-triplet 
spin-1/2 fermions whose left- and right-handed chiral components have the same transformation proper¬ 
ties under the weak-isospin SU(2) gauge group [9, 10]. In these models vector-like quarks are expected to 
couple preferentially to third-generation quarks [9, 1 1] and they can have flavour-changing neutral current 
decays, in addition to the charged-current decays characteristic of chiral quarks. As a result, an up-type 
quark T with charge +2/3 can decay not only to a IT boson and a fi-quark, but also to a Higgs or Z boson 
and a top quark ( T —> Wb. Zt, and Hi). Similarly, a down-type quark B with charge -1/3 can decay to a 
Higgs or Z boson and a fi-quark, in addition to decaying to a IT boson and a top quark ( B —> IT t. Zb , and 
Hb). In order to be consistent with the results from the precision electroweak measurements, a small mass 
splitting between vector-like quarks belonging to the same SU(2) multiplet is required [12], which forbids 
cascade decays such as T —> WB and leaves direct decays into SM particles as the only possibility. Coup¬ 
lings between the vector-like quarks and the first and second quark generations, although not favoured, 
are not excluded [13, 14]. This leads to a rich phenomenology at the LHC, which the experiments are 
investigating. 

Early searches for the pair production of exotic heavy quarks published by the ATLAS and CMS collab¬ 
orations focused on exclusive decay modes assuming a 100% branching ratio. These include searches for 
Tf W + bW-b [15-18], BB -» ZbZb [19-21], and BB W + tW~t [20, 22, 23]. The limits derived 
from these searches cannot easily be applied to other branching ratio values, due to the potentially large 
expected signal contamination from mixed decay modes. A more general search strategy should consider 
simultaneously all three decay modes, providing a more extensive coverage of possible signal contribu¬ 
tions. In absence of an excess, quasi-model-independent limits would be set in the plane defined by the 
branching ratios to two of the decay modes 1 as a function of the heavy-quark mass. The first search that 
considered all three decay modes in the interpretation of results, performed by the ATLAS Collaboration 
using pp collisions at \/v = 7 TeV, primarily targeted the Tt —» W + bW~b process [24]. Using the full 
dataset collected at y/s = 8 TeV, the ATLAS Collaboration has recently published searches for heavy 
quarks decaying to a Z boson and a third-generation quark [25], and searches for heavy quarks decaying 
predominantly to Wt in events with one lepton and jets [26] and in events with two leptons of the same 

1 The branching ratio to the third decay mode is fully determined by the requirement that the sum of branching ratios equals 
unity. 
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charge or three leptons [27]. In the context of vector-like quarks, these searches are used to probe TT and 
BB production, and the three decay modes are considered in the interpretation of the results. The CMS 
Collaboration has published an inclusive search for TT production [28] resulting from the combination 
of several analyses in lepton-plus-jets and multilepton final states at yfs = 8 TeV. This search set 95% 
confidence level (CL) lower limits on the T quark mass ranging between 690 GeV and 780 GeV for all 
possible values of the branching ratios into the three decay modes. 

The results presented in this paper complete the program of searches for pah - production of vector-like 
quarks decaying into third-generation quarks by the ATLAS Collaboration using the pp dataset collected 
at = 8 TeV. Three separate searches are presented, all of them focused on the pair production of 
vector-like quarks in final states involving one isolated electron or muon, high missing transverse mo¬ 
mentum from the undetected neutrino and multiple jets. The first search, referred to as TT — * Wb+X, 
is optimised for TT production with at least one T —> Wb decay, where the resulting W boson acquires 
a high momentum from the large T quark mass. The second search, referred to as TT —> Ht+X, targets 
TT production with at least one T —> Ht decay, with H —> bb, resulting in events with high jet multi¬ 
plicity and a large number of jets tagged as originating from fi-quarks. The third search, referred to as 
BB —> Hb+X, is instead focused on BB production with at least one B —> Hb decay and H —> bb, in 
events with the same final-state signature probed by the TT —> Ht+X search. In all three searches the 
isolated lepton and the high missing transverse momentum are provided by the leptonic decay of a W 
boson originating in the decay of a vector-like quark, a top quark, or a Higgs boson. 

The large mass of the top quark makes it a prime candidate to help uncover the dynamics behind elec- 
troweak symmetry breaking and/or new physics at the electroweak scale. In many new physics models the 
top quark plays a prominent role, often participating in new interactions related to electroweak symmetry 
breaking, or preferentially coupling to new degrees of freedom. Such BSM scenarios usually predict an 
enhanced rate of events containing four top quarks (tttt) in the final state, compared to the SM production 
via the strong interaction. Examples include top quark compositeness [29-31], Randall-Sundrum extra 
dimensions [32], models with coloured scalars [33-38], or universal extra dimensions [39-41]. The CMS 
Collaboration has performed a search for SM tttt production at yfs = 8 TeV in the lepton-plus-jets final 
state [42], setting an observed (expected) 95% CL upper limit on the production cross section of 32 fb 
(32 fb). Using multilepton final states, the ATLAS Collaboration has also searched for SM tttt production 
at = 8 TeV, setting an observed (expected) 95% CL upper limit of 70 fb (27 fb) [27]. The observed 
limit is higher than the expected one owing to an excess of data above the background expectation with 
a significance of 2.5 standard deviations. In addition, the ATLAS multilepton search sensitively probes 
several of the above BSM scenarios giving rise to large enhancements in tttt production. Given its sensit¬ 
ivity to a wide range of models, the TT —» Ht+X search presented in this paper is also used to search for 
a tttt signal, within the SM as well as in the same BSM scenarios as the ATLAS multilepton search, with 
comparable sensitivity. 


2. ATLAS detector 

The ATLAS detector [43] consists of the following main subsystems: an inner tracking system, elec¬ 
tromagnetic and hadronic calorimeters, and a muon spectrometer. The inner detector provides tracking 
information from silicon pixel and microstrip detectors in the pseudorapidity 2 range |//| < 2.5 and from 

2 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis coinciding with the axis of the beam pipe. The .r-axis points from the IP to the centre of the LHC ring, and the 
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a straw-tube transition radiation tracker covering |//| < 2.0, all immersed in a 2 T axial magnetic field 
provided by a superconducting solenoid. The electromagnetic (EM) sampling calorimeter uses lead as 
the absorber material and liquid-argon (LAr) as the active medium, and is divided into barrel (|//| < 1.475) 
and end-cap (1.375 < \rj\ < 3.2) regions. Hadron calorimetry is also based on the sampling technique, 
with either scintillator tiles or LAr as the active medium, and with steel, copper, or tungsten as the ab¬ 
sorber material. The calorimeters cover |//| < 4.9. The muon spectrometer measures the deflection of 
muons with \p\ < 2.7 using multiple layers of high-precision tracking chambers located in a toroidal 
field of approximately 0.5 T and 1 T in the central and end-cap regions of ATLAS, respectively. The 
muon spectrometer is also instrumented with separate trigger chambers covering I 77 I < 2.4. A three-level 
trigger system [44] is used to select interesting events. The first-level trigger is implemented in custom 
electronics and uses a subset of detector information to reduce the event rate to at most 75 kHz. This 
is followed by two software-based trigger levels exploiting the full detector information and yielding a 
typical recorded event rate of 400 Hz during 2012. 


3. Object reconstruction 

The main reconstructed objects considered in this search are electrons, muons, jets, b -jets and missing 
transverse momentum. 

Electron candidates [45] are reconstructed from energy deposits (clusters) in the EM calorimeter that 
are matched to reconstructed tracks in the inner detector. The candidates are required to have a trans¬ 
verse energy 3 Ej greater than 25 GeV and |7/ c iusterl < 2.47, where |? 7 c iusterl is the pseudorapidity of 
the cluster associated with the electron candidate. Candidates in the EM calorimeter transition region 
1.37 < |? 7 ciusterl < 1-52 are excluded. Electrons are required to satisfy “tight” quality requirements [45], 
which include stringent selection requirements on calorimeter, tracking and combined variables that 
provide good separation between prompt electrons and jets. The longitudinal impact parameter of the 
electron track with respect to the event’s primary vertex (see section 4), zo, is required to be less than 2 
mm. To reduce the background from non-prompt electrons resulting from semileptonic decays of b- or 
c-hadrons, and from jets with a high fraction of their energy deposited in the EM calorimeter, electron can¬ 
didates must also satisfy calorimeter- and track-based isolation requirements. The calorimeter isolation 
variable is based on the energy sum of cells within a cone of radius A R = -\/(A 0) * 2 + (A 77) 2 = 0.2 around 
the direction of each electron candidate, and an ^-dependent requirement is made, giving an average effi¬ 
ciency of 90% for prompt electrons from Z boson decays. This energy sum excludes cells associated with 
the electron cluster and is corrected for leakage from the electron cluster itself and for energy deposits 
from additional pp interactions within the same bunch crossing (“pileup”). A further 90%-efficient isol¬ 
ation requirement is made on the track transverse momentum (pj) sum around the electron in a cone of 
radius A R - 0.3. 

Muon candidates [46, 47] are reconstructed from track segments in the various layers of the muon spec¬ 
trometer and matched with tracks found in the inner detector. The final candidates are refitted using 
the complete track information from both detector systems and are required to satisfy p-\ > 25 GeV 

//-axis points upward. Cylindrical coordinates (r,<f>) are used in the transverse plane, (p being the azimuthal angle around the 

beam pipe. The pseudorapidity is defined in terms of the polar angle 9 as 17 = - In tan(0/2). For the purpose of the fiducial 
selection, this is calculated relative to the geometric centre of the detector; otherwise, it is relative to the reconstructed primary 
vertex of each event. 

3 The electron transverse energy is defined as Ej = £ c i us t er / cosh 77 , rack , where Cluster is the energy of the cluster in the calorimeter 
and r/ tr a C k is the pseudorapidity of its associated track. 


5 




and |//| < 2.5. Muons are required to have a hit pattern in the inner detector consistent with a well- 
reconstructed track to ensure good pj resolution. The longitudinal impact parameter of the muon track 
with respect to the primary vertex, zo, is required to be less than 2 mm. Muons are required to satisfy a 
px-dependent track-based isolation requirement: the scalar sum of the pj of the tracks within a cone of 
variable radius A R = 10 GeV//2' around the muon (excluding the muon track itself) must be less than 5% 
of the muon pj (piy). This requirement has good signal efficiency and background rejection even under 
high-pileup conditions, as well as in boosted configurations where the muon is close to a jet. For muons 
from W decays in simulated ft events the average efficiency of the isolation requirement is about 95%. 

Jets are reconstructed with the anti-k, algorithm [48-50] with a radius parameter R = 0.4 from calibrated 
topological clusters [51, 52] built from energy deposits in the calorimeters. Prior to jet finding, a local 
cluster calibration scheme [53] is applied to correct the topological cluster energies for the effects of 
non-compensating response of the calorimeter, dead material and out-of-cluster leakage. The corrections 
are obtained from simulations of charged and neutral particles. After energy calibration [54], jets are 
required to have pj > 25 GeV and |//| < 2.5. To reduce the contamination due to jets originating from 
pileup interactions, a requirement that the so-called “jet vertex fraction" (JVF) be above 0.5 is applied to 
jets with pj < 50 GeV and |?/| < 2.4. This requirement ensures that at least 50% of the scalar sum of 
the pj of the tracks matched to the jet comes from tracks originating from the primary vertex. During 
jet reconstruction, no distinction is made between identified electrons and jet energy deposits. Therefore, 
if any of the jets lie within A R = 0.2 of a selected electron, the closest jet is discarded in order to avoid 
double-counting of electrons as jets. Finally, any electron or muon within A R = 0.4 of a selected jet is 
discarded. 

Jets are identified as originating from the hadronisation of a //-quark (/Magged) via an algorithm [55] that 
uses multivariate techniques to combine information from the impact parameters of displaced tracks as 
well as topological properties of secondary and tertiary decay vertices reconstructed within the jet. For 
each jet, a value for the multivariate /Magging discriminant is calculated, and is considered /Magged if 
this value is above a given threshold. The threshold used in this search corresponds to 70% efficiency 
to tag a //-quark jet, with a light-jet rejection factor 4 of ~ 130 and a charm-jet rejection factor of 5, as 
determined for jets with pj > 20 GeV and |//| < 2.5 in simulated it events. 

The missing transverse momentum (£™ ss ) is constructed [56] from the vector sum of all calorimeter 
energy deposits 5 contained in topological clusters. All topological cluster energies are corrected using 
the local cluster calibration scheme discussed above. Those topological clusters associated with a high- 
pj object (e.g. jet or electron) are further calibrated using their respective energy corrections. In addition, 
contributions from the pj of selected muons are included in the calculation of is™ 88 . 


4. Data sample and event preselection 


This search is based on pp collision data at = 8 TeV collected by the ATLAS experiment between 
April and December 2012. Only events recorded with a single-electron or single-muon trigger under 
stable beam conditions and for which all detector subsystems were operational are considered. The cor¬ 
responding integrated luminosity is 20.3±0.6 fb 1 [57]. Single-lepton triggers with different pj thresholds 

4 The rejection factor is defined as the reciprocal of the selection efficiency. 

5 Each cluster in the calorimeter is considered a massless object and is assigned the four-momentum (// c i us ter, //cluster), where 
//cluster is the measured energy and //cluster is a vector of magnitude //cluster directed from (x, y, z ) = (0,0,0) to the centre of the 
cluster. 
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are combined in a logical OR in order to increase the overall efficiency. The pj thresholds are 24 or 
60 GeV for the electron triggers and 24 or 36 GeV for the muon triggers. The triggers with the lower pj 
threshold include isolation requirements on the candidate lepton, resulting in inefficiencies at high pj that 
are recovered by the triggers with higher pj threshold. Events satisfying the trigger selection are required 
to have at least one reconstructed vertex with at least five associated tracks with pj > 400 MeV, consistent 
with originating from the beam collision region in the x-y plane. The average number of pp interactions 
per bunch crossing is approximately 20, resulting in several vertices reconstructed per event. If more than 
one vertex is found, the hard-scatter primary vertex is taken to be the one which has the largest sum of the 
squared transverse momenta of its associated tracks. For the event topologies considered in this paper, 
this requirement leads to a probability to reconstruct and select the correct hard-scatter primary vertex 
larger than 99%. 

Events are required to have exactly one reconstructed electron or muon and at least four jets satisfying the 
quality and kinematic criteria discussed in section 3. The selected lepton is required to match, with A R < 
0.15, the lepton reconstructed by the trigger. The background from multijet production is suppressed 
by a requirement on E™ 1SS as well as on the transverse mass of the lepton and £™ ss (/??” )/’ For both 
lepton selections the requirements are £'" 11SS > 20 GeV and £ , ' | 11ISS + nij > 60 GeV. Further suppression 
of the background not including ft-quark jets is achieved by requiring at least one ^-tagged jet in the 
TT —» Wb+X search, and at least two ^-tagged jets in the TT —> Ht+X and BB —> Hb+X searches. 
In the following, events satisfying either the electron or muon selections arc combined and treated as a 
single analysis channel. 


5. Signal modelling 


This section describes the different signal scenarios considered in the interpretation of the results, together 
with details of how they are modelled in the analysis. 


5.1. Vector-like quark pair production 

Vector-like quarks with mass below approximately 1 TeV are mostly produced in pairs via the strong inter¬ 
action in pp collisions at yfs = 8 TeV. For higher masses, single production mediated by the electroweak 
interaction can potentially dominate, depending on the strength of the interaction between the new quarks 
and the weak gauge bosons. The predicted pair-production cross section ranges from 5.3 pb for a quark 
mass of 350 GeV to 3.3 fb for a quark mass of 1000 GeV, with an uncertainty that increases from 8% to 
14% over this mass range. This cross section is independent of the electroweak quantum numbers of the 
new heavy quark and just depends on its mass. It was computed using Top++ v2.0 [58] at next-to-next-to- 
leading order (NNFO) in QCD, including resummation of next-to-next-to-leading logarithmic (NNFF) 
soft gluon terms [59-63], and using the MSTW 2008 NNFO [64, 65] set of parton distribution functions 
(PDF). Theoretical uncertainties result from variations on the factorisation and renormalisation scales, as 
well as from uncertainties on the PDF and as- The latter two represent the largest contribution to the over¬ 
all theoretical uncertainty on the cross section and were calculated using the PDF4FHC prescription [66] 
with the MSTW 2008 68% CL NNLO, CT10 NNLO [67, 68] and NNPDF2.3 5f FFN [69] PDF sets. 


6 mjr = 


yflp'- 


cos A (j>), where p r T is the transverse momentum (energy) of the muon (electron) and A0 is the azimuthal 
angle separation between the lepton and the direction of the missing transverse momentum. 


t ,E™ ss (1 
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As discussed in section 1, vector-like quarks can couple preferentially to third-generation quarks, as 
the mixing between weak eigenstates of the same electric charge is proportional to the mass of the SM 
quark [9, 11], and thus present a rich phenomenology. In particular, a vector-like quark has neutral- 
current tree-level decays to a Z or H boson plus a SM quark, in addition to the charged-current decay 
mode to a W boson and a SM quark, which is the only decay mode chiral quarks can have. Figure 1 
depicts representative Feynman diagrams for the signals probed by the searches discussed in this paper. 
The branching ratios to each of these decay modes vary as a function of the heavy-quark mass and depend 
on its weak-isospin (SU(2)) quantum numbers [10]. Figure 2(a) shows the branching ratios as a function 
of T quark mass in the SU(2) singlet and doublet hypotheses. 7 In the case of a singlet, all three decay 
modes have sizeable branching ratios, while the charged-current decay mode T —» Wb is absent in the 
doublet cases. The doublet prediction is valid for an (X, T ) doublet, where the charge of the X quark is 
+5/3, as well as a ( T , B) doublet when a mixing assumption of \Vjb\ ^ IU/;I is made, where V\j are the 
elements of a generalised Cabibbo-Kobayashi-Maskawa matrix [10]. Since the T quark branching ratios 
are identical in both doublets, in the following no distinction between them is made when referring to the 
T quark doublet hypothesis. Similarly, figure 2(b) shows the branching ratios as a function of B quark 
mass in the singlet and doublet hypotheses. In the case of a ( T , B) doublet with the mixing assumption 
\Vrb\ ^ I V,b|, BR(B —» Wt) = 1, while such a decay mode is absent for the {B, Y ) doublet case, where the 
charge of the Y quark is -4/3. The Y quark is equivalent to a chiral quark since it only has charged-current 
decays, Y —» W~b. 

Simulated samples of TT and BB are generated with the leading-order (LO) generator Protos v2.2 [70] 
using the MSTW 2008 LO PDF set and passed to Pythia 6.426 [7 1] for parton showering and fragmenta¬ 
tion. The AUET2B [72] set of optimised parameters for the underlying event (UE) description, referred to 
as the “UE tune”, is used. The vector-like quarks are forced to decay with a branching ratio of 1 /3 to each 
of the three modes (W, Z, H). Arbitrary sets of branching ratios consistent with the three decay modes 
summing to unity are obtained by reweighting the samples using particle-level information. Samples are 
generated assuming singlet couplings and for heavy-quark masses between 350 GeV and 1100 GeV in 
steps of 50 GeV. Additional samples are produced at two mass points (350 GeV and 600 GeV) assuming 
doublet couplings in order to confirm that kinematic differences arising from the different chirality of 
singlet and doublet couplings are negligible in this analysis. In all simulated samples (both signal and 
background) used in this search, the top quark and SM Higgs boson masses are set to 172.5 GeV and 
125 GeV respectively. The samples are normalised using the Top++ cross section predictions discussed 
above. 


5.2. Four-top-quark production 

The production cross section for four-top-quark events in the SM is very small (ergg — I fb at yfs = 
8 TeV) [73, 74], but it can be significantly enhanced in several BSM scenarios. Figure 3 depicts rep¬ 
resentative LO Feynman diagrams for four-top-quark production within the SM and the different BSM 
scenarios considered in this paper. A class of models involving new heavy vector particles strongly 
coupled to the right-handed top quark, such as top quark compositeness [29-31] or Randall-Sundrum ex¬ 
tra dimensions [32], can be described via an effective field theory (EFT) involving a four-fermion contact 


7 The branching ratios in figure 2 are valid for small mixing between the new heavy quark and the third-generation quark [ 10- 
12 ]. 





Figure 1: Representative leading-order Feynman diagrams for TT production probed by (a) the TT —> Wb+X 
search and (b) the TT —> Ht+X search, and (c) for BB production probed by the BB —> Hb+X search. 




Figure 2: Branching ratios for the different decay modes as a function of heavy-quark mass in the case of (a) a 
vector-like T quark and (b) a vector-like B quark, as computed with Protos. In both cases the branching ratios are 
provided for an SU(2) singlet and two different SU(2) doublet scenarios. 
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Figure 3: Representative leading-order Feynman diagrams for four-top-quark production within (a) the SM and 
several beyond-the-SM scenarios (see text for details): (b) via an effective four-top-quark interaction in an effective 
field theory model, (c) via scalar-gluon-pair production, and (d) via cascade decays from Kaluza-Klein excitations 
in a universal extra dimensions model with two extra dimensions compactified using the geometry of the real 
projective plane. 


interaction [75] (figure 3(b)). The Lagrangian assumed is 

■C-At = \-^(tRy IJ tR)(tR'y t jtR), (1) 

where /r is the right-handed top quark spinor, y^ are the Dirac matrices, C 4 , is the coupling constant, and 
A is the energy scale of new physics. Only the contact interaction operator with right-handed top quarks 
is considered, since left-handed operators are already strongly constrained by the precision electroweak 
measurements [76], 

In addition, two specific models are considered involving new heavy particles: scalar gluon (sgluon) 
pair production and a Universal Extra Dimensions (UED) model. Sgluons are colour-adjoint scalars, 
denoted by cr, that appear in several extensions of the SM, both supersymmetric [33, 34] and non- 
supersymmetric [35-38]. The dominant production mode at the LHC is in pairs via the strong interaction, 
gg —> crcr. For sgluon masses above twice the top quark mass, the dominant decay mode is into tt, giving 
rise to a four-top-quark final state (figure 3(c)). The UED model considered has two extra dimensions that 
are compactified using the geometry of the real projective plane (2UED/RPP) [39], leading to a discret¬ 
isation of the momenta along their directions. A tier of Kaluza-Klein towers is labelled by two integers, 
k and €, referred to as “tier (k, €)”. Within a given tier, the squared masses of the particles are given at 
leading order by m 2 = k 2 /R 2 + ( 2 /R 2 , where nR+ and nR*, are the size of the two extra dimensions. The 
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model is parameterised by R 4 and R$ or, alternatively, by //?kk = I /R 4 and = A 4 /R 5 . Four-top-quark 
production can arise from tier ( 1 , 1 ), where particles from this tier have to be pair produced because of 
symmetries of the model. Then they chain-decay to the lightest particle of this tier, the heavy photon 
A (1,1) , by emitting SM particles (figure 3(d)). The branching ratios of into SM particles are not 
predicted by the model, although the decay into tt is expected to be dominant [40]. Four-top-quark events 
can also arise from tiers (2,0) and (0,2) via a similar mechanism. In this case the expected cross section 
for four-top-quark production is reduced compared to that from tier ( 1 , 1 ) since each state in tiers ( 2 , 0 ) and 
(0,2) can decay directly into a pair of SM particles or into a pair of states in tiers (1,0) or (0,1) via bulk in¬ 
teractions, resulting in smaller branching ratios for decay into tt [40]. In the following, when considering 
four-top-quark production from a given tier, it is assumed that the A photon in that tier decays with 100% 
branching ratio into tt while A photons from other tiers cannot decay into tt. Within this model, observa¬ 
tions of dark-matter relic abundance prefer values of »7 kk between 600 GeV and 1200 GeV [41]. 

Simulated samples of four-top-quark production within the SM, within an EFT model, and within the 
2UED/RPP model, are generated with the Madgraph5 1.3.33 [77] LO generator and the MSTW 2008 
PDF set, interfaced to Pythia 8.1 [78] and the AU2 UE tune [79]. In the case of the 2UED/RPP model, 
samples are generated for four different values of tmkk (600, 800, 1000 and 1200 GeV) and the Bridge [80] 
generator is used to decay the pair-produced excitations from tier (1,1) generated by Madgraph5. Con¬ 
straints for tiers ( 2 , 0 ) and ( 0 , 2 ) can be derived from those for tier ( 1 , 1 ) together with the theoretical cross 
sections. Samples of four-top-quark production via sgluon pairs are generated with Pythia 6.426 with the 
CTEQ6L1 [81] PDF set and the AUET2B UE tune, for seven different values of the sgluon mass between 
350 GeV and 1250 GeV, and normalised to the NLO theoretical cross section [82], 

Events from minimum-bias interactions are simulated with the Pythia 8.1 generator with the MSTW 
2008 LO PDF set and the A2 tune [79]. They are overlaid on the simulated signal events according to 
the luminosity profile of the recorded data. The contributions from these pileup interactions are modelled 
both within the same bunch crossing as the hard-scattering process and in neighbouring bunch crossings. 
Finally, the generated samples arc processed through a simulation [83] of the detector geometry and 
response using Geant4 [84] with a fast simulation of the calorimeter response [83]. Ah samples arc 
processed through the same reconstruction software as the data. Simulated events are corrected so that 
the object identification efficiencies, energy scales and energy resolutions match those determined from 
data control samples. 


6. Background modelling 


After event preselection, the main background is ff+jets production, with the production of a W boson 
in association with jets (H'+jcts) and multijet events contributing to a lesser extent. Small contributions 
arise from single top quark, Z+jets and diboson (WW, WZ , ZZ) production, as well as from the associated 
production of a vector boson V (V = W, Z) or a Fliggs boson and a tt pair (ttV and ttH). Multijet events 
contribute to the selected sample via the misidentification of a jet or a photon as an electron or via the 
presence of a non-prompt lepton, e.g. from a semileptonic b- or c-hadron decay; the corresponding yield 
is estimated via data-driven methods. The rest of the background contributions are estimated from simula¬ 
tion and normalised to then - theoretical cross sections. In the case of the ff+jets and W/Z +jets background 
predictions, further corrections are applied to improve agreement between the data and simulation, as 
discussed in sections 6.1 and 6.2 respectively. 
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All simulated background samples utilise Photos 2.15 [85] to simulate photon radiation and Tauola 
1.20 [86] to simulate r decays. Similarly to the signal samples, they also include a simulation of pileup 
interactions, and are processed through a full Geant4 detector simulation and the same reconstruction 
software as the data. Further details about the modelling of each of the backgrounds are provided be¬ 
low. 


6.1. G+jets background 

Simulated samples of ff+jets events are generated with the next-to-leading-order (NLO) generator Powheg- 
Box 2.0 [87-90] using the CT10 PDF set [67]. The nominal sample is interfaced to Pythia 6.425 [71] with 
the CTEQ6L1 PDF set and the Perugia2011C UE tune [91]. An alternative sample, used to study the un¬ 
certainty related to the fragmentation model, is interfaced to FIerwig v6.520 [92] with the CTEQ6L1 PDF 
set and Jimmy v4.31 [93] to simulate the UE. The ff+jets samples are normalised to the theoretical cross 
section obtained with Top++, performed at NNLO in QCD and including resummation of NNLL soft 
gluon terms. 

The ff+jets samples are generated inclusively, but events are categorised depending on the flavour content 
of additional particle jets in the event (i.e. jets not originating from the decay of the tt system). Particle 
jets are reconstructed with the anti-k, algorithm with a radius parameter R = 0.4 and arc required to have 
pj > 15 GeV and |//| < 2.5. Events where at least one such particle jet is matched within A R < 0.4 to a 
/i-hadron with pj > 5 GeV not originating from a top quark decay are labelled as tt + bb events. Similarly, 
events where at least one such particle jet is matched within A R < 0.4 to a c-hadron with pj > 5 GeV not 
originating from a W boson decay, that are not labelled already as tt + bb, are labelled as tt + cc events. 
Events labelled as either tt + bb or tt + cc are generically referred to below as ff+HF events, where F1F 
stands for “heavy flavour”. The remaining events are labelled as ff+light-jet events, including those with 
no additional jets. In Powheg+Pythia the modelling of ff+HF is via the parton-shower evolution. To 
study uncertainties related to this simplified description, an alternative fl+jets sample is generated with 
Madgraph5 1.5.11 using the CT10 PDF set. It includes tree-level diagrams with up to three additional 
partons (including b- and c-quarks) and is interfaced to Pythia 6.425. 

Since the best possible modelling of the fl+jets background is a key aspect of these searches, a correction 
is applied to simulated tt events in Powheg+Pythia based on the ratio of the differential cross sections 
measured in data and simulation at ^fs = 7 TeV as a function of top quark pj and tt system p-\ [94]. The 
stability of the ratio between yfs = 1 TeV and yfs = & TeV was studied to support the usage of yfs = 
7 TeV data to correct the simulation at = 8 TeV. This correction significantly improves agreement 
between simulation and data in distributions such as the jet multiplicity and the pj of decay products 
of the tt system. This correction is applied only to ff+light-jets and tt + cc events. The modelling of 
the tt + bb background, particularly important for the Ht/Hb+X searches, is improved by reweighting 
the Powheg+Pythia prediction to an NLO prediction of tt + bb including parton showering [95], based 
on Sherpa+OpenLoops [96, 97] using the CT10 PDF set. This reweighting is performed for different 
topologies of tt + bb in such a way that the inter-normalisation of each of the categories and the relevant 
kinematic distributions are at NLO accuracy. More details about the modelling of the ff+jets background 
can be found in ref. [98], 
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6.2. W/Z +jets background 


Samples of W/Z+jets events are generated with up to five additional partons using the Alpgen v2. 14 [99] 
LO generator and the CTEQ6L1 PDF set, interfaced to Pythia v6.426 for parton showering and fragment¬ 
ation. To avoid double-counting of partonic configurations generated by both the matrix-element calcu¬ 
lation and the parton shower, a parton-jet matching scheme (“MLM matching”) [100] is employed. The 
VV'+jcts samples are generated separately for W'-t-light-jets, Wbb+ jets, VTcc+jets, and VVA+jets. The Z+jets 
samples are generated separately for Z+light-jets, Zbb+ jets, and Zee +jets. Overlap between VQQ +jets 
(V = W,Z and Q = b,c) events generated from the matrix-element calculation and those generated from 
parton-shower evolution in the VV'/Z+light-jets samples is avoided via an algorithm based on the angular 
separation between the extra heavy quarks: if A R(Q, Q) > 0.4, the matrix-element prediction is used, 
otherwise the parton-shower prediction is used. Both the VV'+jcts and Z+jets background contributions 
are normalised to their inclusive NNLO theoretical cross sections [101]. Further corrections are applied 
to W/Z+jets events in order to better describe data in the preselected sample. Scale factors for each 
of the VV'+jcts categories (VV7A+jcts, VTcc+jets, VVA+jcts and W'+l ight-jets) are derived for events with 
one lepton and at least four jets by simultaneously analysing six different event categories, defined by 
the b-tag multiplicity (0, 1 and >2) and the sign of the lepton charge. The b-tag multiplicity provides 
information about the heavy-flavour composition of the W +jets background, while the lepton charge is 
used to determine the normalisation of each component, exploiting the expected charge asymmetry for 
VV'+jcts production in pp collisions as predicted by Alpgen. In the case of Z+jets events, a correction to 
the heavy-flavour fraction was derived to reproduce the relative rates of Z+2-jets events with zero and 
one ^-tagged jets observed in data. In addition, the Z boson pj spectrum was compared between data 
and the simulation in Z+2-jets events, and a reweighting function was derived in order to improve the 
modelling. 


6.3. Other simulated background 

Samples of single-top-quark backgrounds corresponding to the /-channel, ^-channel and Wt production 
mechanisms are generated with Powheg-Box 2.0 [102, 103] using the CT10 PDF set and interfaced to 
Pythia 6.425 with the CTEQ6F1 PDF set and the Perugia2011C UE tune. Overlaps between the tt and 
Wt final states are removed using the “diagram removal” scheme [104]. The single-top-quark samples are 
normalised to the approximate NNFO theoretical cross sections [105-107] calculated using the MSTW 
2008 NNFO PDF set. 

The WW/WZ/ZZ +jets samples are generated with up to three additional partons using Alpgen v2.13 and 
the CTEQ6F1 PDF set, interfaced to Herwig v6.520 and Jimmy v4.3 1 for parton showering, fragmentation 
and UE modelling. The MEM parton-jet matching scheme is used. The VV'VV'+jcts samples require at 
least one of the W bosons to decay leptonically, while the VV'Z/ZZ+jcts samples require one Z boson to 
decay leptonically, with the other boson decaying inclusively. Additionally, WZ+jets samples requiring 
the W and Z bosons to decay leptonically and hadronically respectively, are generated with up to three 
additional partons (including massive b- and c-quarks) using Sherpa vl.4.1 and the CT10 PDF set. All 
diboson samples are normalised to then - NFO theoretical cross sections [108] 

Samples of ttV events, including ttWW, are generated with up to two additional partons using Madgraph5 
1.3.28 with the CTEQ6F1 PDF set, and interfaced to Pythia 6.425 with the AUET2B UE tune. A sample 
of tlH events is generated with the PowHel framework [109], which combines the Powheg-Box generator 
and NFO matrix elements obtained from the HEFAC-Oneloop package [110]. The sample is generated 
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using the CTlOnlo PDF set [67]. Showering is performed with Pythia 8.1 using the CTEQ6L1 PDF set 
and the AU2 UE tune [72, 111]. Inclusive decays of the Higgs boson are assumed in the generation of the 
ttH sample. The ttV samples are normalised to the NLO cross section predictions [1 12]. The ttH sample 
is normalised using the NLO cross section [113-115] and the Higgs decay branching ratios [116-119] 
collected in ref. [120]. 


6.4. Multijet background 

Multijet events can enter the selected data sample through several production and misreconstruction 
mechanisms. In the electron channel, the multijet background consists of non-prompt electrons as well as 
misidentified photons (e.g. with a conversion into an e + e~ pair) or jets with a high fraction of their energy 
deposited in the EM calorimeter. In the muon channel, the background contributed by multijet events is 
predominantly due to final states with non-prompt muons, such as those from semileptonic b- or c-hadron 
decays. 

The multijet background normalisation and shape are estimated directly from data by using the “matrix 
method” technique [121]. The matrix method exploits differences in lepton-identification-related proper¬ 
ties between prompt, isolated leptons from W and Z boson decays (referred to as “real leptons” below) 
and those where the leptons are either non-isolated or result from the misidentification of photons or jets 
(referred to as “fake leptons” below). For this puipose, two samples are defined after imposing the final 
kinematic selection criteria, differing only in the lepton identification criteria: a “tight” sample and a 
“loose” sample, the former being a subset of the latter. The tight selection employs the complete set of 
lepton identification criteria used in the analysis. For the loose selection the lepton isolation requirements 
are omitted. The method assumes that the number of selected events in each sample (/V loose and /V llghl ) 
can be expressed as a linear combination of the numbers of events with real and fake leptons, so that the 
number of multijet events in the tight sample is given by 

<f = £fake (frealN l0 °‘ Se - N^), (2) 

treat £fake 

where e rea ] (tfake) represents the probability for a real (fake) lepton that satisfies the loose criteria to 
also satisfy the tight ones. Both probabilities are measured in data control samples. To measure e rea |, 
samples enriched in real leptons from W bosons decays are selected by requiring high /:' | n ' ss or m^. 
The average e rea t is -0.75 (-0.98) in the electron (muon) channel. To measure samples enriched in 
multijet background are selected by requiring either low £™ ss (electron channel) or high impact parameter 
significance for the lepton track (muon channel). The average epake value is -0.35 (-0.20) in the electron 
(muon) channel. Dependencies of e rea | and ejakc on quantities such as lepton pj and //, A R between the 
lepton and the closest jet, or number of ^-tagged jets, are parameterised in order to obtain a more accurate 
estimate. 

7. Search for TT —> Wb+X production 

This search is sensitive to TT production where at least one of the T quarks decays into a W boson and a 
b-quark, although it is particularly optimised for TT —> W + bW~b events. One of the W bosons present in 
the final state is then required to decay leptonically. After the preselection described in section 4, further 
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background suppression is achieved by applying requirements aimed at exploiting the distinct kinematic 
features of the signal. The large T quark mass results in energetic W bosons and ^-quarks in the final 
state with large angular separation between them, while the decay products from the boosted W bosons 
have small angular separation. The combination of these properties is very effective in distinguishing the 
dominant tt background since tt events with boosted W boson configurations are rare and are typically 
characterised by a small angular separation between the W boson and the 6 -quark from the top quark 
decay. 

To take advantage of these properties, it is necessary to identify the hadronically decaying W boson (Whad) 
as well as the b- jets in the event. The candidate /;-jcts are defined as the two jets with the highest 6 -tag 
discriminant value, although only one of them is explicitly required to be 6 -tagged in the event selec¬ 
tion. Two types of Wf ia d candidates are defined, W^ y d e 1 and 11 , depending on the angular separation 
between their decay products. W^ y d e 1 candidates correspond to boosted W bosons, where the quarks from 
the W-boson decay emerge with small angular separation and are reconstructed as a single jet. Alternat¬ 
ively, 11 candidates are characterised by two reconstructed jets. In the construction of both types of 
IThad candidates, the two candidate 6-jets are not considered. 

A W d ype 1 candidate is defined as a single jet with pj > 400 GeV, which is the typical pj above which 
the decay products from a W boson would have an angular separation A R < tfcone = 0.4. A < P d eH 
candidate is defined as a dijet system with pj > 250 GeV, angular separation A R(j, j) < 0.8 and mass 
within the range of 60-120 GeV. The asymmetric window about the VV'-boson mass value is chosen in 
order to increase the acceptance for hadronically decaying Z bosons from TT —* WbZt events. Any jets 
satisfying the W* ype 1 requirements are excluded from consideration when forming W ^ d e11 candidates. 
The leptonically decaying W boson (Wj e p) is reconstructed using the lepton and £™ ss , which is taken 
as a measurement of the neutrino pj. Requiring that the invariant mass of the lepton-neutrino system 
equals the nominal W boson mass allows reconstruction of the longitudinal momentum of the neutrino 
up to a two-fold ambiguity. If two solutions exist, they are both considered. If no real solution exists, the 
pseudorapidity of the neutrino is set equal to that of the lepton, since in the kinematic regime of interest 
the decay products of the W boson tend to be collinear. 

Table 1 summarises the event selection requirements. Two selections, “loose” and “tight”, are defined, 
with the latter being more restrictive than the former and representing the final selection. As discussed 
below, the loose selection is used to validate the background modelling in a kinematic regime close to 
the final selection. The loose selection considers preselected events with at least one W ^ 61 or W^ y d e 11 
candidate. If multiple Whad candidates are found in a given event, the one with the highest p\ is chosen. 
Figure 4(a) shows the distribution of the number of Whad candidates after preselection. The events must 
satisfy Hj > 800 GeV, where Hj is the scalar sum of the lepton pj, T n [ uss and the pj of the selected 
jets. The Hj distribution peaks at ~2mj for signal events, which makes the Hj > 800 GeV requirement 
particularly efficient for signal with mj > 400 GeV, while rejecting a large fraction of the background. 
Figure 4(b) shows the distribution of Hj after the requirement of >1 Whad candidate and prior to the 
Hj > 800 GeV requirement. In addition, the highest-px 6 -jet candidate (b\ ) and the next-to-highest-px 
6 -jet candidate (hi) are required to have pj(b\) > 160 GeV and pj{bT) > 110(80) GeV respectively, 
in the case of a W* y d eI (W' y d eH ) candidate. Finally, the angular separation between the lepton and the 
reconstructed neutrino is required to satisfy A R{£, v) < 0.8 ( 1 . 2 ) in case of a W^ y d e 1 (W d y d e n ) candidate. 
Figure 5(a) shows the distributions of A R(£,v) after all previous requirements and prior to the AR(C, v) 
requirement. 

The tight selection adds further requirements that are particularly effective at suppressing tt background. 
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Selection 

Requirements 

Preselection 

Exactly one electron or muon 
£™ iss > 20 GeV, £™ ss + > 60 GeV 

>4 jets, >1 fi-taggedjets 

Loose selection 

Preselection 

> 1 IFhad candidate (type I or type II) 

Hj > 800 GeV 

Pt(^i) > 160 GeV, pj{bi) > 110 GeV (type I) or pxfe) > 80 GeV (type II) 

A R{b, v) < 0.8 (type I) or A R(t, v) < 1.2 (type II) 

Tight selection 

Loose selection 

min(AR(£, b\p)) > 1.4, min(A/?(IFhad> ^ 1 , 2 )) > 1-4 

AR(b\,b 2 ) > 1.0 (type I) or AR{b\,b 2 ) > 0.8 (type II) 

Am < 250 GeV (type I) [see text for definition] 


Table 1: Summary of event selection requirements for the TT —> Wb+X analysis (see text for details). 


First, a large angular separation between the W bosons and the b -jets from the top quark decay is re¬ 
quired: min(A/?(£,£ 1 , 2 )) > 1.4 and min(A/?(IFhad>Zfi, 2 )) > 1.4. Figure 5(b) shows the distributions of 
min(Ai?(7’, ^ 12 )) after loose selection and prior to the min(A/?(^, b\p)) >1.4 requirement. Finally, addi¬ 
tional requirements are made on AR{b\,b 2 ) > 1.0(0.8) in the case of a IF^ 0 ) candidate and 

Am < 250 GeV only in the case of a 1 candidate, where Am = mindm'^o - /n^J) is the smallest 
absolute difference between the reconstructed heavy-quark masses obtained by pairing the fFi ep and IFhad 
candidates with the two b -jet candidates as described in the following. The reconstruction of the VF| ep 
candidate usually yields two solutions, and there are two possible ways to pair - the Z?-jet candidates with 
the W boson candidates to form the heavy quarks. Among all possible combinations, the one yielding the 
smallest Am is chosen. The main discriminating valuable used in this search is the reconstructed heavy- 
quark mass (m ieco ), built from the IFhad candidate and one of the two b-jct candidates. The resulting m rec0 
distributions for the loose and tight selections are shown in figure 6 for the sum of IF^ e1 and W^ e 11 
events. The tight selection has the better expected sensitivity, and only this selection is chosen to derive 
the final result of the search. The loose selection, displaying a significant ff background at low m reco which 
is in good agreement with the expectation, provides further confidence in the background modelling prior 
to the application of fi-jet isolation requirements in the tight selection. 

Table 2 presents a summary of the background estimates for the loose and tight selections, as well as a 
comparison of the total predicted and observed yields. The quoted uncertainties include both the statistical 
and systematic contributions. The latter are discussed in section 10. The predicted and observed yields 
are in agreement within these uncertainties. 
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Figure 4: TT —* Wb+X search: distribution of (a) the number of hadronically decaying W boson (Wh a d) candidates 
after preselection requirements, and (b) the scalar sum ( Hj ) of the transverse momenta of the lepton, the selected jets 
and the missing transverse momentum after preselection and >1 Whad candidate requirements. The data (solid black 
points) are compared to the SM prediction (stacked histograms). The contributions from backgrounds other than 
tt are combined into a single background source referred to as “Non-ff”. The total uncertainty on the background 
estimation is shown as a black hashed band. The expected contribution from a vector-like T quark with mass 
mj = 600 GeV under the assumption BR(T —> Wb) = 1, multiplied by a factor of 50, is also shown (red dashed 
histogram). The lower panel shows the ratio of data to the SM prediction. The last bin contains the overflow. 


8. Search for TT —> Ht+X and tttt production 


This search is focused on Tf production where at least one of the T quarks decays into a Higgs boson 
and a top quark resulting from the following processes: TT —* HtHt, ZtHt and WbHt: For the dominant 
H —* bb decay mode, the final-state signature is characterised by high jet and b-tag multiplicities, which 
provide a powerful experimental handle to suppress the background. Similarly, this search is also sensitive 
to Tf —> ZtZt and WbZt, with Z —> bb. High jet and h-tag multiplicities are also characteristic of tttt 
events (both within the SM and in BSM extensions), which makes this search also sensitive to this process. 
Figure 7(a) compares the jet multiplicity distribution after preselection (described in section 4) between 
the total background and several signal scenarios. Signal events have, on average, higher jet multiplicity 
than the background. The higher ft-quark content of signal events results in a higher Z?-tag multiplicity 
than for the background, as illustrated in figure 7(b) for events with >6 jets. Therefore, after preselection, 
the final selection requirements are >5 jets of which >2 jets are ^-tagged, leaving a sample completely 
dominated by fH-jets background. In order to ensure a non-overlapping analysis sample and to facilitate 
the combination of results, events accepted by the Wb+X search are rejected. This veto only removes 


In the following ZtHt is used to denote both ZtHt sad its charge conjugate, HtZt. Similar notation is used for other processes, 
as appropriate. 
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AR(Z,v) min(AR(Z,b i 2 )) 
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Figure 5: TT —> 147;+X search: distribution of (a) the angular separation between the lepton and the reconstruc¬ 
ted neutrino (A R{(, v)), and (b) the minimum angular separation between the lepton and the two candidate /?-jets 
(min(A/?(£, ^ 1 , 2 )))- The selections made include all previous requirements except for the requirement on each of 
these variables (see text for details). The data (solid black points) are compared to the SM prediction (stacked his¬ 
tograms). The contributions from backgrounds other than tt are combined into a single background source referred 
to as “Non-ff”. The total uncertainty on the background estimation is shown as a black hashed band. The expected 
contribution from a vector-like T quark with mass mj = 600 GeV under the assumption BR(T —> Wh) - 1 is also 
shown (red histogram), stacked on top of the SM background. The lower panel shows the ratio of data to the SM 
prediction. The last bin contains the overflow. 


about 2% of the events with >6 jets and >4 ^-tagged jets in data. 

In order to optimise the sensitivity of the search, the selected events are categorised into different channels 
depending on the number of jets (5 and >6) and on the number of ^-tagged jets (2, 3 and >4). The 
channel with >6 jets and >4 ^-tagged jets has the largest signal-to-background ratio and therefore drives 
the sensitivity of the search. The channels with two and three ^-tagged jets have significantly lower 
signal-to-background ratio. These are particularly useful to calibrate the ft+jets background prediction 
and constrain the related systematic uncertainties. In the case of the channel with >6 jets and >4 ^-tagged 
jets the background uncertainty is dominated by uncertainties on the ^-tagging, jet energy calibration and 
physics modelling, including the ff+HF content. A detailed discussion of the systematic uncertainties 
considered is given in section 10. In addition, events with >6 jets and 3 or >4 ^-tagged jets are split 
into two channels each depending on the value of the invariant mass of the two ^-tagged jets with lowest 
A R separation: M™ nAR < 100 GeV (“low M™ nA *”) and M™ nAS > 100 GeV (“high M™ inAS ”). For high 
values of mj, the Higgs boson from the T —> Ht decay has high pj, and the bb pair from the Higgs 
boson decay has smaller angular separation than other pairs resulting from combinatorial background. As 
shown in figure 8(a), the variable provides a good approximation to the reconstructed H —* bb 

invariant mass and allows the separation of these channels into channels depleted or enriched in T —» Ht, 
H —» bb decays, the latter having a higher signal-to-background ratio. Therefore, the total number of 
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Figure 6: TT —» Wb+X search: distribution of the reconstructed heavy-quark mass (m reco ) after (a) the loose selec¬ 
tion and (b) the tight selection, for the sum of 1 and W ^ 11 events. The data (solid black points) are compared 
to the SM prediction (stacked histograms). The contributions from backgrounds other than tt are combined into 
a single background source referred to as “Non-tt”. The total uncertainty on the background estimation is shown 
as a black hashed band. The expected contributions from a vector-like T quark with mass m r = 600 GeV in two 
scenarios, BR(T —> Wb) = 1 (red histogram) and singlet (dashed black histogram), are also shown stacked on top 
of the SM background. The lower panel shows the ratio of data to the SM prediction. The last bin contains the 
overflow. 
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Figure 7: TT —> Ht+X search (simulated events): comparison of (a) the jet multiplicity distribution after preselec¬ 
tion, and (b) the b -tag multiplicity distribution after the requirement of >6 jets, between the total background 
(shaded histogram) and several signal scenarios considered in this search: TT production in the T quark singlet 
(red solid histogram) and doublet (red dashed histogram) cases, and sgluon pair production giving a four-top-quark 
final state (red dotted histogram). A mass of 600 GeV is assumed for the T quark and the sgluon. 
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Tight selection 

TT (m T = 600 GeV) 



BRfT ->Wb)=l 

115 + 10 

58.9 + 5.9 

Singlet 

60.3 + 5.1 

24.5 + 2.3 


it 

390+ 110 

10.7 + 4.3 

tiV 

6.5 ±2.5 

0.4 ± 0.2 

ttH 

1.6 ±0.4 

0.10 + 0.03 

VF+jcts 

38+19 

11.4 + 6.2 

Z+jets 

1.5 + 1.2 

0.4 ± 0.4 

Single top 

36+17 

2.2+ 1.5 

Diboson 

5.6+ 1.4 

1.5 ±0.6 

Multijet 

0.3+ 1.6 

0.8 ±0.7 


Total background 

480+ 120 

27.6 + 8.6 

Data 

478 

34 


Table 2: TT —» Wh+X search: number of observed events, integrated over the whole mass spectrum, compared to 
the SM expectation after the loose and tight selections. The expected signal yields in two different scenarios for a 
vector-like T quark with m T = 600 GeV, BR(T —> Wb) = 1 and singlet, are also shown. The quoted uncertainties 
include both the statistical and systematic contributions. 


analysis channels considered in this search is eight: (5 j, 2 b), (5 j, 3 b), (5 j, >4 b), (>6 j, 2 b), (>6 j, 3 b, 
low M™ nAS ), (>6 j, 3 b, high M™ nAS ), (>6 j, >4 b, low M™ nAS ), and (>6 j, >4 b, high M™ nAR ), where 
(n j, m b) indicates n selected jets and m /^-tagged jets. 

To further improve the separation between signal and background, the distinct kinematic features of the 
signal are exploited. In particular, the large T quark mass results in energetic leptons and jets in the final 
state, and Hj provides a suitable discriminating variable between signal and background. Figure 8(b) 
compares the Hj distribution between signal and background for events with >6 jets and >4 fi-tagged 
jets. The Hj distribution is quite similar for different signal scenarios corresponding to pair production 
of exotic particles with the same mass (600 GeV in this case), and significantly different from that of the 
background. The discrimination between signal and background increases with mass. 

Figures 9 and 10 show the comparison of data and prediction for the Hj distributions in each of the ana¬ 
lysis channels considered. The corresponding predicted and observed yields per channel can be found in 
table 3. Following the statistical procedure outlined in section 11, a fit to the observed Hj distributions 
in data in the eight analysis channels is performed. This provides an improved background prediction 
with smaller uncertainties, and hence improved sensitivity to a signal. The results are presented in sec¬ 
tion 12. 


9. Search for BB —> Hb+X production 

This search is focused on BB production where at least one of the B quarks decays into a Higgs boson and 
a b quark, a decay mode that was omitted from previous searches [25-27], In particular, the BB —» HbHb 
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5j,2b 

5j, 3 b 

5 j, >4 b 

>6 j, 2 b 

TT (m T = 600 GeV) 





Singlet 

52.5 ±4.2 

19.0 ±2.3 

5.8 ± 1.2 

123.3 ± 6.2 

(T, B ) or (X, T) doublet 

25.8 ±2.0 

14.0 ± 1.4 

5.0 ± 1.0 

154.1 ±6.4 

era —> tttt (mu- = 800 GeV) 

2.0 ±0.3 

1.4 ±0.3 

0.3 ±0.1 

64.8 ±4.6 

tttt+X (Tier (1,1), «?kk = 800 GeV) 

1.0 ±0.4 

0.6 ±0.3 

0.06 ± 0.05 

180 ±29 

ff+light-jets 

32400 ± 5300 

2930 ± 520 

48 ± 12 

16200 ± 4000 

tt + cc 

3800 ±2100 

730 ±410 

42 ±24 

3300 ± 1800 

tt + bb 

1530 ± 800 

800 ± 420 

108 ±58 

1300 ± 700 

ttV 

140 ± 46 

24.9 ± 8.1 

2.9 ± 1.0 

172 ± 56 

ttH 

39.2 ± 1.7 

20.8 ± 1.6 

5.6 ±0.7 

60.2 ± 4.5 

W+jets 

1600± 1000 

111 ±71 

5.0 ±3.4 

770 ± 530 

Z+jets 

360± 120 

24.8 ± 8.4 

1.2 ±0.5 

185 ±67 

Single top 

1630 ± 320 

169 ± 36 

7.0 ± 1.0 

730 ± 200 

Diboson 

85 ±27 

7.3 ±2.5 

0.4 ±0.2 

45 ± 15 

Multijet 

133 ±48 

33 ± 12 

6.9 ± 2.6 

56 ±20 

Total background 

41700 ±6400 

4840 ± 900 

228 ± 69 

22800 ± 5200 

Data 

43319 

5309 

244 

23001 



>6 j, 3 b 
low M“ inAS 

>6 j, 3 b 
high M™ nAS 

>6 j, >4 b 
low M“ inAS 

>6 j, >4 b 
high M™ aAR 

TT (m T = 600 GeV) 





Singlet 

29.5 ± 2.0 

44.0 ± 3.6 

17.7 ± 1.9 

24.1 ±3.7 

(T, B) or (V, T) doublet 

50.2 ±2.5 

68.9 ±4.1 

41.0 ±3.9 

53.8 ±7.3 

crcr —* tttt (»v = 800 GeV) 

22.5 ± 1.6 

50.7 ± 3.5 

9.3 ± 1.0 

16.2 ± 2.6 

tttt+X (Tier (1,1), = 800 GeV) 

33.6 ±2.8 

132.5 ±5.9 

27.7 ± 2.3 

75 ± 13 

/f+light-jets 

1280 ±350 

440± 110 

38 ± 14 

9.3 ±3.9 

tt + cc 

550 ± 320 

220 ± 120 

53 ±31 

14.7 ±9.0 

tt + bb 

620 ± 330 

250± 140 

178 ±95 

46 ±25 

ttV 

28.7 ±9.2 

12.5 ±4.2 

6.2 ± 2.0 

1.5 ±0.5 

ttH 

24.9 ± 1.9 

11.6 ± 1.3 

10.6 ± 1.2 

4.1 ±0.6 

lV±jets 

68 ±46 

16 ± 10 

6.6 ±4.8 

0.6 ±0.4 

Z±jets 

15.7 ±6.3 

3.3 ± 1.3 

1.6 ± 0.6 

0.3 ±0.1 

Single top 

74 ±22 

32 ± 12 

7.8 ±2.2 

2.1 ± 1.3 

Diboson 

4.2 ± 1.6 

1.2 ±0.5 

0.4 ±0.1 

0.2 ± 0.1 

Multijet 

1.9 ±0.8 

4.8 ±2.1 

< 0.01 

2.8 ± 1.0 

Total background 

2670 ± 680 

990 ± 260 

300 ± 110 

81 ±30 

Data 

3015 

1085 

362 

84 


Table 3: TT —* Ht+X search: predicted and observed yields in each of the analysis channels considered. The 
background prediction is shown before the fit to data. Also shown are the signal predictions for different benchmark 
scenarios considered. The quoted uncertainties are the sum in quadrature of statistical and systematic uncertainties 
on the yields. 
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Figure 8: TT —» Ht+X search (simulated events): comparison of the distributions of (a) the invariant mass of the two 
^-tagged jets with lowest A R separation and (b) the scalar sum of the transverse momenta of the lepton, 

the selected jets and the missing transverse momentum ( Hj ), between the total background (shaded histogram) and 
several signal scenarios considered in this search: TT —> WbHt (red solid histogram), TT —> WbZt or SM tttt 
production (red dashed histograms), and sgluon pair production giving a tttt final state (red dotted histogram). A 
mass of 600 GeV is assumed for the T quark and the sgluon. The selection used in both (a) and (b) corresponds to 
events satisfying the preselection requirements and with >6 jets and >4 ^-tagged jets. 


final state is the least covered one because the most-common Higgs boson decay mode, H —> bb, leads 
to a challenging final state with six b- jets and no leptons. In contrast, cleaner experimental signatures 
involving leptons tend to be suppressed by the small decay branching ratios. However, a sizeable signal 
rate results from the mixed decay mode where one of the Higgs bosons decays into W + W~, while the 
other Higgs boson decays into bb\ BB —» HbHb —> (W + W~)b(bb)b. When one of the W bosons decays 
leptonically, this leads to the final-state signature considered in this search, involving one lepton and high 
jet and b- tag multiplicities, analogous to the signature exploited by the TT —» Ht+X search. 

Consequently, this search considers the same discriminating variable, Hj, and the same eight analysis 
channels as the TT —» Ht+X search. Figure 11(a) illustrates the good separation between signal and 
background in the Hj distribution for events passing the preselection requirements and with > 6 jets and 
> 4 ^-tagged jets. A peculiarity of the B —> Hb decay mode is that the fi-jct originating (directly) from 
the B-quark decay can have very high transverse momentum in the case of a heavy B quark. To exploit 
this feature, the event selection is tightened relative to that used in the TT —> Ht+X search by raising the 
minimum pj requirement on the two highest-pT (leading) ^-tagged jets to pj > 150 GeV. Figure 11(b) 
shows the distribution of the subleading b- jet pj for events passing the preselection requirements and 
with > 6 jets and > 4 ^-tagged jets. The tighter requirement on the subleading A-jct pj rejects about 
90% of the tt background while retaining a large acceptance for the BB —> Hb+X signal. This search 
is also sensitive to other BB final states, such as BB —> HbWt, that typically do not involve multilepton 
final states in the topologies usually searched for (opposite-sign dileptons with aZ -> ( + £~ candidate, 
same-sign dileptons, and trileptons), and is thus complementary to previous searches [25-27], 

Figures 12 and 13 show the comparison of data and prediction for the H\ distributions in each of the 
analysis channels considered. The corresponding predicted and observed yields per channel can be found 
in table 4. The results of the fit to the data to improve the background prediction, as in the TT —> Ht+X 
search, are presented in section 12. 
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Figure 9: TT —> Ht+X search: comparison between data and prediction for the distribution of the scalar sum 
(7/i ) of the transverse momenta of the lepton, the selected jets and the missing transverse momentum in each of 
the analysed channels after final selection: (a) (5 j, 2 b), (b) (5 j, 3 b), (c) (5 j, >4 b), and (d) (>6 j, 2 b). The 
background prediction is shown before the fit to data. The contributions from W/Z+jets, single top, diboson and 
multijet backgrounds are combined into a single background source referred to as “Non-ff”. Also shown is the 
expected signal contribution from a singlet vector-like T quark with mass m T - 600 GeV. The last bin in all figures 
contains the overflow. The bottom panel displays the ratio of data to the total background prediction. The hashed 
area represents the total uncertainty on the background. 
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Figure 10: TT —> Ht+X search: comparison between data and prediction for the distribution of the scalar sum 
(Hj) of the transverse momenta of the lepton, the selected jets and the missing transverse momentum in each of 
the analysed channels after final selection: (a) (>6 j, 3 b, low M££ nAR ), (b) (>6 j, 3 b, high M^ nAR ), (c) (>6 j, >4 
b, low and (d) (>6 j, >4 b, high M^ nAR ). The background prediction is shown before the fit to data. The 

contributions from W/Z+ jets, single top, diboson and multijet backgrounds are combined into a single background 
source referred to as “Non-ff”. Also shown is the expected signal contribution from a singlet vector-like T quark 
with mass hit = 600 GeV. The last bin in all figures contains the overflow. The bottom panel displays the ratio of 
data to the total background prediction. The hashed area represents the total uncertainty on the background. 
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5 j, 2 b 

5 j, 3 b 

5 j, >4 b 

>6 j, 2 b 

BB (m B = 600 GeV) 





BR(B ->Hb) = 1 

8 .6+ 1.1 

9.3 ±2.2 

5.0+ 1.4 

11.9 + 3.0 

Singlet 

12.2+ 1.9 

8.8 ± 1.7 

3.4 ±0.8 

27.4 + 4.3 

(B, F) doublet 

8.5 ± 1.1 

5.8 + 1.4 

2.8 ±0.8 

10.9 + 2.1 

/f+light-jets 

389 ± 93 

72+ 18 

2.1 ±0.7 

234 ± 74 

tt + cc 

56 + 42 

23 + 15 

2.2+ 1.5 

55 + 40 

tt + bb 

19+ 14 

25 ± 14 

5.5 ±3.2 

22+ 15 

ttV 

4.2+ 1.4 

1.6 ±0.5 

0.3 ±0.1 

5.1 ± 1.7 

ttH 

1.0 ±0.1 

1.1 ±0.2 

0.5 ±0.1 

1.5 ±0.2 

IF+jets 

21 ± 12 

3.5 ±2.1 

0.6 ±0.5 

12.5 + 7.9 

Z+jets 

8.2 ±3.3 

2.8 ±2.8 

0.5 ±0.5 

4.3 ±4.1 

Single top 

41.3 + 7.2 

8.8 + 1.9 

0.6 ±0.1 

28.0 + 6.8 

Diboson 

1.9 ±0.9 

0.5 ±0.3 

0.07 ± 0.07 

1.2 ±0.7 

Multijet 

<0.01 

<0.01 

0.4 ± 0.2 

0.2 ±0.1 

Total background 

540 ± 120 

139 ± 35 

12.8 + 4.9 

360 ± 100 

Data 

576 

165 

10 

375 



>6 j, 3 b 
low M“ inAS 

>6 j, 3 b 
high M™ nAR 

>6 j, >4 b 
low M“ inAS 

>6 j, >4 b 
high M“ inAS 

BB (m B = 600 GeV) 
BR(B -* Hb) = 1 
Singlet 
(B, F) doublet 

3.8 ±0.6 

7.1 ±0.9 

2.7 ±0.3 

13.1 ± 1.8 
15.8 + 2.5 

7.0+ 1.3 

3.2 ±0.7 
4.6 ±0.9 

2.3 ±0.6 

9.6 ±2.0 

7.5 + 1.5 

3.9 ±0.9 

ff+light-jets 

21.3 + 9.0 

32.8 + 9.5 

1.4 ±0.5 

1.5 ±0.6 

tt + cc 

10.8 + 7.5 

20+15 

2 .2+ 1.6 

2.9 ±2.2 

tt + bb 

13.1 + 8.5 

24+ 16 

7.8 ±4.8 

8.1 ±5.3 

ttV 

1.1 ±0.4 

1.6 ±0.6 

0.6 ± 0.2 

0.4 ± 0.2 

ttH 

0.7 ±0.1 

1.4 ±0.2 

0.5 ±0.1 

0.9 ±0.2 

VF+jets 

2.0+ 1.3 

1.1 ±0.8 

0.3 ±0.3 

0.05 ± 0.05 

Z+jets 

0.11+0.07 

0.2 ±0.1 

<0.01 

<0.01 

Single top 

3.2 ±0.6 

5.1 ±2.2 

0.8 ±0.2 

0.3 ±0.2 

Diboson 

0.2 ±0.1 

0.09 ± 0.03 

0.02 ± 0.01 

<0.01 

Multijet 

< 0.01 

0.6 ± 0.2 

<0.01 

0.4 ±0.1 

Total background 

53 + 18 

87 + 30 

13.7 + 5.9 

14.5 + 7.3 

Data 

62 

103 

23 

20 


Table 4: BB —> Hb+X search: predicted and observed yields in each of the analysis channels considered. The 
background prediction is shown before the fit to data. Also shown are the signal predictions for different benchmark 
scenarios considered. The quoted uncertainties are the sum in quadrature of statistical and systematic uncertainties 
on the yields. 


25 



> 

CD 

0 

O 

o 


0.4- 
0.35 F 
0.31- 
0.25^ 

0.2 j- 

0.15 ^ 
0.1 H 
0.05^ 


ATLAS 

Simulation 
Vs = 8 TeV 


JL 


| | Total background 

- BB -> HbHb (600 GeV) 

. BB singlet (600 GeV) 

. BB (B,Y) doublet (600 GeV) 

> 6 jets, > 4 b-tags 


500 




1000 


1500 2000 

H t [GeV] 



(a) 
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Figure 11: BB —> Hb+X search (simulated events): comparison of the distributions of (a) the scalar sum of the trans¬ 
verse momenta of the lepton, the selected jets and the missing transverse momentum (H\), and (b) the transverse 
momentum of the next-to-highest-transverse-momentum /;-jet, between the total background (shaded histogram) 
and several BB signal scenarios considered in this search: BR(B —» lib) = 1 (red solid histogram), B quark singlet 
(red dashed histogram), and B quark from a ( B , Y) doublet (red dotted histogram). In all cases a mass of 600 GeV 
is assumed for the B quark. The selection used in both (a) and (b) corresponds to events satisfying the preselection 
requirements and with >6 jets and >4 ft-tagged jets. 
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(C) (d) 

Figure 12: BB —> Hb+X search: comparison between data and prediction for the distribution of the scalar sum 
(7/ r ) of the transverse momenta of the lepton, the selected jets and the missing transverse momentum in each of the 
analysed channels after final selection: a) (5 j, 2 b), (b) (5 j, 3 b), (c) (5 j, >4 b), and (d) (>6 j, 2 b). The background 
prediction is shown before the fit to data. The contributions from W/Z+ jets, single top, diboson and multijet 
backgrounds are combined into a single background source referred to as “Non-ff”. Also shown is the expected 
signal contribution from a vector-like B quark with mass ms = 600 GeV under the assumption BRffi —» I lb) = 1. 
The last bin in all figures contains the overflow. The bottom panel displays the ratio of data to the total background 
prediction. The hashed area represents the total uncertainty on the background. 
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Figure 13: BB —> Hb+X search: comparison between data and prediction for the distribution of the scalar sum 
( Hj ) of the transverse momenta of the lepton, the selected jets and the missing transverse momentum in each of 
the analysed channels after final selection: (a) (>6 j, 3 b, low M^ mAS ), (b) (>6 j, 3 b, high (c) (>6 j, >4 

b, low M'™ nAK ), and (d) (>6 j, >4 b, high M^ nAR ). The background prediction is shown before the fit to data. The 
contributions from W/Z+ jets, single top, diboson and multijet backgrounds are combined into a single background 
source referred to as “Non-ff”. Also shown is the expected signal contribution from a vector-like B quark with 
mass ms = 600 GeV under the assumption BR( B —> I lb) = 1. The last bin in all figures contains the overflow. 
The bottom panel displays the ratio of data to the total background prediction. The hashed area represents the total 
uncertainty on the background. 
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10. Systematic uncertainties 


Several sources of systematic uncertainty are considered that can affect the normalisation of signal and 
background and/or the shape of their corresponding final discriminant distributions. Individual sources 
of systematic uncertainty are considered uncorrelated. Correlations of a given systematic uncertainty 
are maintained across processes and channels. Table 5 presents a list of all systematic uncertainties 
considered in the analyses and indicates whether they arc taken to be normalisation-only, or to affect both 
shape and normalisation. 

Table 6 presents a summary of the systematic uncertainties for the TT —> Wb+X search and then - impact 
on the normalisation of signal and backgrounds. A similar summary is presented for the TT —* Ht+X 
and BB —> Hb+X searches in tables 7 and 8 respectively, restricted to the highest-sensitivity channel and 
displaying only the signal and the tt+jets background categories. Tables 7 and 8 also show the impact of 
the systematic uncertainties before and after the fit to data. 

In the case of the TT —> Wb+X search, the total systematic uncertainty in the background normalisation 
is approximately 29%, with the dominant contributions originating from the normalisation of the W'+jcts 
background (20%), jet energy scale (+17%/-12%) and the tt+HF normalisation (11%). The total system¬ 
atic uncertainty in the signal normalisation is +8%/-10%, with comparable contributions from jet energy 
scale and ^-tagging uncertainties. 

The leading sources of systematic uncertainty in the TT —> Ht+X and BB —> Hb+X searches vary de¬ 
pending on the analysis channel considered, but they typically originate from tt +jets modelling (including 
tf+HF), jet energy scale and ^-tagging. For example, the total systematic uncertainty in the background 
normalisation in the highest-sensitivity channel (>6 j, >4 b, high M^ nAR ) of the TT —> Ht+X search 
is approximately 37%, with the largest contributions originating from f?+HF normalisation (23%), jet 
energy scale (10%) and ^-tagging (9%). However, as discussed previously, the fit to data in the eight 
analysis channels in these searches allows the overall background uncertainty to be reduced significantly, 
to approximately 5% in the case of the TT —> Ht+X search. More details about the fit to data can be 
found in section 12.1. The total systematic uncertainty on the signal normalisation is approximately 15%, 
almost all due to ^-tagging uncertainties. 

The following sections describe each of the systematic uncertainties considered in the analyses. 


10.1. Luminosity 

The uncertainty on the integrated luminosity is 2.8%, affecting the overall normalisation of all processes 
estimated from the simulation. It is derived following the same methodology as that detailed in ref. [57], 

10.2. Reconstructed objects 
10.2.1. Leptons 

Uncertainties associated with leptons arise from the reconstruction, identification and trigger, as well as 
the lepton momentum scale and resolution. The reconstruction and identification efficiency of electrons 
and muons, as well as the efficiency of the trigger used to record the events, differ slightly between data 
and simulation. Scale factors are derived using tag-and-probe techniques onZ -> £ + (t = e,jT) data 
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Systematic uncertainty 

Type 

Components 

Luminosity 

N 

1 

Reconstructed Objects 



Electron 

SN 

5 

Muon 

SN 

6 

Jet reconstruction 

SN 

1 

Jet vertex fraction 

SN 

1 

Jet energy scale 

SN 

22 

Jet energy resolution 

SN 

1 

Missing transverse momentum 

SN 

2 

^-tagging efficiency 

SN 

6 

c-tagging efficiency 

SN 

4 

Light-jet tagging efficiency 

SN 

12 

High-px tagging 

SN 

1 

Background Model 



tt cross section 

N 

1 

tt modelling: p-\ re weighting 

SN 

9 

tt modelling: parton shower 

SN 

3 

ff+HF: normalisation 

N 

2 

tt+cc : HF re weighting 

SN 

2 

tt+cc: generator 

SN 

4 

tt+bb: NFO Shape 

SN 

8 

IL+jcts normalisation 

N 

3 

Z+jets normalisation 

N 

3 

Single top cross section 

N 

1 

Single top model 

SN 

1 

Diboson normalisation 

N 

1 

ttV cross section 

N 

1 

ttV model 

SN 

1 

ttH cross section 

N 

1 

ttH model 

SN 

2 

Multijet normalisation 

N 

2 


Table 5: List of systematic uncertainties considered. An “N" means that the uncertainty is taken as normalisation- 
only for all processes and channels affected, whereas “SN" means that the uncertainty is taken on both shape and 
normalisation. Some of the systematic uncertainties are split into several components for a more accurate treatment. 


and simulated samples, and are applied to the simulation to correct for differences. Additional sources of 
uncertainty originate from the corrections applied to adjust the lepton momentum scale and resolution in 
the simulation to match those in data, measured using reconstructed distributions of the Z —> t + f ~ and 
J/i/j —» t + (~ masses, as well as the measured E/p in W —> ev events, where E and p arc the electron 
energy and momentum, as measured by the calorimeter and the tracker respectively. The combined effect 
of all these uncertainties results in an overall normalisation uncertainty on the signal and background of 
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Signal 

tt 

Non-tf 

Total background 

Luminosity 

±2.8 

±2.8 

±2.8 

±2.8 

Lepton efficiencies 

±1.6 

±1.6 

±1.5 

±1.6 

Jet energy scale 

±3.4/-7.2 

±16 

+ 19/-9 

±17/-12 

Jet efficiencies 

±1.5 

±1.6 

±1.6 

±1.6 

Jet energy resolution 

±1.1 

±0.6 

±2.6 

±1.8 

6-tagging efficiency 

±5.0 

±0.7 

±2.9 

±2.0 

c-tagging efficiency 

±0.4 

±1.2 

±2.3 

±1.9 

Light-jet tagging efficiency 

±0.2 

±1.3 

±1.6 

±1.4 

High-/;- r tagging efficiency 

±3.2 

±1.3 

±0.8 

±1.1 

Missing transverse momentum 

- 

±2.6 

- 

±1.0 

tt : re weighting 

- 

±15 

- 

±5.9 

tv. parton shower 

- 

±9.3 

- 

±3.6 

tf+HF: normalisation 

- 

+12.0/—5.5 

- 

+4.5/—2.1 

tf+HF: modelling 

- 

±30 

- 

±11 

Theoretical cross sections 

- 

±6.0 

±33 

±20 

Multijet normalisation 

- 

- 

±2.9 

±1.8 

Non-tf modelling 

- 

- 

±2.3 

±1.4 

Total 

+7.7/-10.0 

±40 

±35 

±29 


Table 6: TT —> Wb+X search: summary of the systematic uncertainties considered and their impact (in %) on 
the normalisation of signal and backgrounds. Only sources of systematic uncertainty resulting in a normalisation 
change of at least 0.5% are displayed. The signal shown corresponds to a vector-like T quark with mass m T = 
600 GeV and BR(T -> Wb) = 1. 


approximately 1.5%. 


10.2.2. Jets and missing transverse momentum 

Uncertainties associated with jets arise from the efficiency of jet reconstruction and identification based 
on the JVF variable, as well as the jet energy scale and resolution. The uncertainty associated with the jet 
reconstruction efficiency is assessed by randomly removing 0.2% of the jets with pj below 30 GeV, which 
is the level of disagreement between data and the simulation, and has a negligible impact in the analysis. 
The per-jet efficiency to satisfy the JVF requirement is measured in Z(—> ( + t ~)+1 -jet events in data 
and simulation, selecting separately events enriched in hard-scatter jets and events enriched in jets from 
pileup, and good agreement is found. The associated uncertainty is estimated by changing the nominal 
JVF cut value by ±0.1 and repeating the analysis using the modified cut value, resulting in normalisation 
uncertainties in the range of 1-5%, depending on the jet multiplicity under consideration and the p-\ 
spectra of the jets. The jet energy scale and its uncertainty were derived by combining information from 
test-beam data, LHC collision data and simulation [54]. The jet energy scale uncertainty is split into 22 
uncorrelated sources with their respective jet pj and // dependences and are treated independently in this 
analysis. It represents one of the leading sources of uncertainty associated with reconstructed objects, 
affecting the normalisations of signal and backgrounds by approximately 5% and 15% respectively, in 
the most signal-rich channels considered. The jet energy resolution was measured in data and simulation 
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>6 j, >4 b, high 



Signal 

Pre-fit 

tH-light-jets 

tt + cc 

tt + 66 

Post-fit 

tH-light-jets tt + cc 

tt + 66 

Luminosity 

± 2.8 

± 2.8 

± 2.8 

± 2.8 

± 2.6 

± 2.6 

± 2.6 

Lepton efficiencies 

±1.5 

±1.5 

±1.5 

±1.5 

±1.5 

±1.4 

±1.5 

Jet energy scale 

±4.4 

±15 

±11 

±12 

±8.7 

±6.4 

±6.7 

Jet efficiencies 

- 

±4.0 

± 2.2 

±1.9 

±2.7 

±1.5 

±1.3 

Jet energy resolution 

± 0.1 

±4.4 

±3.8 

±0.5 

±3.1 

± 2.6 

±0.4 

6 -tagging efficiency 

±13 

±5.6 

±5.4 

±9.3 

±4.6 

±4.6 

± 6.6 

c-tagging efficiency 

± 1.6 

±5.8 

±12 

±3.1 

±5.6 

±11 

±2.9 

Light-jet tagging efficiency 

± 0.6 

±20 

±5.7 

± 2.0 

±17 

±5.1 

± 1.8 

High-pT tagging efficiency 

±4.8 

±0.7 

±1.7 

± 1.6 

± 0.6 

±1.3 

± 1.2 

tt: reweighting 

- 

±13 

±15 

- 

±10 

±10 

- 

tt: parton shower 

- 

±28 

±17 

± 6.2 

±13 

±11 

±4.0 

77+HL: normalisation 

- 

- 

±50 

±50 

- 

±32 

±18 

ff±HL: modelling 

- 

- 

±17 

±12 

- 

±16 

±10 

Theoretical cross sections 

- 

±6.3 

±6.3 

±6.3 

±4.6 

±4.6 

±4.6 

Total 

±15 

±42 

±61 

±55 

±22 

±30 

±15 


Table 7: TT —» Ht+X search: summary of the systematic uncertainties considered in the (>6 j, >4 b, high M]™ nAS ) 
channel and their impact (in %) on the normalisation of signal and backgrounds, before and after the fit to data. 
Only sources of systematic uncertainty resulting in a normalisation change of at least 0.5% are displayed. The 
signal shown corresponds to a singlet vector-like T quark with mass m T = 600 GeV. The total post-fit uncertainty 
can be different from the sum in quadrature of individual sources due to the anti-correlations between them resulting 
from the fit to the data. 


as a function of jet pj and rapidity using dijet events. They are found to agree within 10%, and the 
corresponding uncertainty is assessed by smearing the jet pj in the simulation. 

The £'™ ss reconstruction is affected by uncertainties associated with leptons and jet energy scales and 
resolutions, which are propagated to fi™ ss and thus are included under the corresponding uncertainty 
categories in tables 6 - 8 . Additional small uncertainties associated with the modelling of the underlying 
event, in particular its impact on the pj scale and resolution of unclustered energy, are also taken into 
account and are displayed in tables 6-8 under the category of “Missing transverse momentum". 


10.2.3. Heavy- and light-flavour tagging 

Efficiencies to tag jets from 6 - and c-quarks in the simulation are corrected to match the efficiencies 
in data by py-dependent factors in the approximate ranges 0.9-1.0 and 0.9-1.1 respectively, whereas 
the light-jet efficiency is scaled by pj- and 77 -dependent scale factors in the range 1.2-1.5 [55, 122]. 
Uncertainties on these scale factors include a total of six independent sources affecting 6 -jets and four 
independent sources affecting c-jets. Each of these uncertainties has different jet pj dependence. Twelve 
uncertainties are considered for the light-jets tagging, which depend on the jet pj and 77 regions. These 
systematic uncertainties are taken as uncorrelated between 6 -jets, c-jets, and light-jets. An additional 
uncertainty is included due to the extrapolation of the 6 -, c-, and light-jet-tagging scale factors for jets 
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>6 j, >4 b, high 



Signal 

Pre-fit 

tt+light-jets tt + cc 

tt + 66 

Post-fit 

tt+light-jets tt + cc 

tt + 66 

Luminosity 

+2.8 

±2.8 

±2.8 

±2.8 

±2.7 

±2.7 

±2.7 

Lepton efficiencies 

±1.6 

±1.4 

±1.5 

±1.7 

±1.4 

±1.5 

±1.6 

Jet energy scale 

±5.6 

±14 

±14 

±11 

±13 

±14 

±11 

Jet efficiencies 

±3.1 

±3.3 

±1.0 

±0.9 

±3.2 

±0.9 

±0.8 

Jet energy resolution 

±0.1 

±6.0 

±1.1 

±1.9 

±4.5 

±0.9 

±1.5 

6 -tagging efficiency 

±16 

±7.6 

±9.2 

±16 

±3.9 

±5.2 

±7.5 

c-tagging efficiency 

±1.0 

±6.1 

±15 

±3.0 

±5.8 

±14 

±2.8 

Light-jet tagging efficiency 

- 

±19 

±6.3 

±2.4 

±18 

±5.8 

±2.3 

High-pT tagging efficiency 

±11 

±2.7 

±5.3 

±5.0 

±1.9 

±3.8 

±3.6 

tt: reweighting 

- 

±15 

±16 

- 

±14 

±15 

- 

tt: parton shower 

- 

±22 

±35 

±26 

±14 

±33 

±24 

tt+HF: normalisation 

- 

- 

±50 

±50 

- 

±44 

±30 

//+HF: modelling 

- 

- 

±27 

±24 

- 

±28 

±21 

Theoretical cross sections 

- 

±6.3 

±6.2 

±6.3 

±5.9 

±5.9 

±5.9 

Total 

±21 

±38 

±73 

±65 

±24 

±46 

±27 


Table 8: BB —» Hb+X search: summary of the systematic uncertainties considered in the (>6 j, >4 b, high M™ nhR ) 
channel and their impact (in %) on the normalisation of signal and backgrounds, before and after the fit to data. 
Only sources of systematic uncertainty resulting in a normalisation change of at least 0.5% are displayed. The 
signal shown corresponds to a vector-like B quark with mass ms = 600 GeV and BR(7i —> lib) - 1. The total 
post-fit uncertainty can be different from the sum in quadrature of individual sources due to the anti-correlations 
between them resulting from the fit to the data. 


with pj beyond the kinematic reach of the data calibration samples used: pj > 300 GeV for b- and c-jets, 
and pj > 750 GeV for light-jets. This uncertainty is evaluated in the simulation by comparing the tagging 
efficiencies while varying e.g. the fraction of tracks with shared hits in the silicon detectors or the fraction 
of fake tracks resulting from random combinations of hits, both of which typically increase at high pj 
due to growing track multiplicity and density of hits within the jet. These uncertainties are taken to be 
correlated among the three jet flavours. As an example, the uncertainties on the tagging efficiencies for 
6 -jets and c-jets with 300 GeV < pj < 500 GeV arc 14% and 23% respectively. 


10.3. Background modelling 
10.3.1. ff+jets 

A number of systematic uncertainties affecting the modelling of tf+jets are considered. These include the 
uncertainty on the theoretical prediction for the inclusive cross section, uncertainties associated with the 
reweighting procedure applied to ?7+light-jets and tt + cc processes, uncertainties affecting the modelling 
of ff+HF-jets production, and uncertainties associated with the choice of parton shower and hadronisa- 
tion model. A summary of these uncertainties can be found below. Additional details can be found in 
ref. [98], 
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An uncertainty of +5%/-6% is assumed for the inclusive tt production cross section [58], including 
contributions from varying the factorisation and renormalisation scales and uncertainties arising from the 
PDF, as and the top quark mass. The PDF and as uncertainties were calculated using the PDF4LF1C 
prescription. 

Uncertainties associated with the reweighting procedure applied to ff+light-jets and tt + cc processes 
include the nine leading sources of uncertainty in the dilferential cross section measurement at \fs = 
7 TeV [94], dominated by the modelling of initial- and final-state radiation and the choice of event gener¬ 
ator for tt production. 

Uncertainties affecting the modelling of tt + bb production include those associated with the NLO pre¬ 
diction from Sherpa+OpenLoops, which is used for reweighting of the default Powheg tt + bb prediction. 
These include three different scale variations, including changing the functional form of the renormal¬ 
isation scale, changing the functional form of the factorisation and resummation scales, and varying the 
renormalisation scale by a factor of two up and down. In addition, a different shower recoil model scheme 
and two alternative PDF sets (MSTW and NNPDF) are considered. A fraction of the tt + bb background 
predicted by Powheg+Pythia originates from multiple parton interactions or final-state radiation from 
top decay products. Such backgrounds are not part of the NLO prediction, and these two categories are 
kept separate and subject to additional normalisation uncertainties. The NLO corrections and associated 
systematic uncertainties are adjusted so that the overall normalisation of the tt + bb background at the 
particle level is fixed, i.e. effectively only migrations across categories and distortions to the shape of 
the kinematic distributions are considered. Detailed comparisons of tt + bb between Powheg+Pythia and 
Sherpa+OpenLoops show that the cross sections agree to better than 50%, which is taken as a normalisa¬ 
tion uncertainty for tt + bb. 

Beyond the uncertainties associated with the reweighting procedure, additional uncertainties are assigned 
to the modelling of the tt + cc component of the background, which again is not part of the NLO pre¬ 
diction used for tt + bb. These include two uncertainties taken as the full difference between applying 
and not applying the reweightings of the top quark and tt p-\ spectra. In addition, four uncertainties arc 
considered associated with the choice of LO generator: the full difference between Powheg+Pythia and 
Madgraph5+Pythia simulations, as well as variations in generator parameters (factorisation and renor¬ 
malisation scales, matching threshold and c-quark mass), which are derived using Madgraph5+Pythia 
simulations and applied to the Powheg+Pythia simulation. Analogously to the procedure used in the 
tt + bb background estimate, these uncertainties are adjusted so that the overall normalisation of the tt + cc 
background at the particle level is fixed. Finally, an overall normalisation uncertainty of 50% is also 
assigned to the tt + cc component, taken as uncorrelated with the same normalisation uncertainty applied 
to tt + bb, since only the tt + bb process is normalised to a NLO prediction. 

An uncertainty due to the choice of parton shower and hadronisation model is derived by comparing 
events produced by Powheg interfaced to Pythia or Herwig. In the case of ff+light-jets and tt + cc, a 
reweighting of the top quark and tt p-\ spectra is also applied to the Powheg+Herwig samples to ensure 
reliable modelling of the top quark kinematics. The corresponding correction factors were recalculated 
for Powheg+Herwig in order to match the differential cross section measurements at yfs = 1 TeV. In the 
case of tt + bb, the various HF categories and the corresponding partonic kinematics in Powheg+Herwig 
are reweighted to match the NLO prediction of Sherpa+OpenLoops, so that only the effect of changing 
the hadronisation model is propagated. Given the different effect of this uncertainty on the ff+light-jets, 
tt + cc and tt + bb, it is treated as uncorrelated between the three processes. This treatment prevents an 
undue reduction of this systematic uncertainty on tt + cc and tt + bb by constraining it for tt+ light-jets via 
the fit to data in the highly populated channels with two fi-tagged jets. 
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10.3.2. W/Z+ jets 


Uncertainties affecting the modelling of the VT/Z+jets background include 5% from their respective nor¬ 
malisations to the theoretical NNLO cross sections [101], as well as an additional 24% normalisation 
uncertainty added in quadrature for each additional inclusive parton multiplicity bin, based on a compar¬ 
ison among different algorithms for merging LO matrix elements and parton showers [123]. The above 
uncertainties are taken as uncorrelated between VV'+jets and Z+jets. 


10.3.3. Other simulated background 

Uncertainties affecting the modelling of the single-top-quark background include a +5%/-4% uncertainty 
on the total cross section estimated as a weighted average of the theoretical uncertainties on Wt- and s- 
channel production [105-107], as well as a systematic uncertainty on FA-channel production concerning 
the separation between tt and Wt at NLO [124]. The latter is estimated by comparing the nominal sample, 
which uses the so-called “diagram subtraction” scheme, with an alternative sample using the “diagram 
removal” scheme. 

Uncertainties on the diboson background normalisation include 5% from the NLO theoretical cross sec¬ 
tions [108] added in quadrature to an uncertainty of 24% due to the extrapolation to the high jet multipli¬ 
city channels, following the procedure discussed in section 10.3.2. 

Uncertainties on the ttV and ttH normalisations are 30% and +9%/-12% respectively, from the uncer¬ 
tainties on their respective NLO theoretical cross sections [1 12, 120, 125]. Additional small uncertainties 
arising from scale variations, which change the amount of initial-state radiation and thus the event kin¬ 
ematics, are also included. 


10.3.4. Multijet 

Uncertainties on the data-driven multijet background estimate receive contributions from the limited 
sample size in data, particularly at high jet and b -tag multiplicities, as well as from the uncertainty on 
the rate of fake leptons, estimated in different control regions (e.g. selected with either an upper £™ 1SS 
or requirement). A combined normalisation uncertainty of 50% due to all these effects is assigned, 
which is taken as correlated across jet and h-tag multiplicity bins, but uncorrelated between electron and 
muon channels. No explicit shape uncertainty is assigned since the large statistical uncertainties associ¬ 
ated with the multijet background prediction, which are uncorrelated bin-to-bin in the final discriminating 
variable, effectively cover all possible shape uncertainties. 


11. Statistical analysis 


For a given search, the distributions of the final discriminating variables in each of the analysis channels 
considered are combined to test for the presence of a signal. The statistical analysis is based on a binned 
likelihood function £,(ja, 6) constructed as a product of Poisson probability terms over all bins considered 
in the analysis. This function depends on the signal-strength parameter /j., a multiplicative factor to the 
theoretical signal production cross section, and 6, a set of nuisance parameters that encode the effect of 
systematic uncertainties on the signal and background expectations and are implemented in the likelihood 
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function as Gaussian or log-normal priors. Therefore, the total number of expected events in a given bin 
depends on p and 6. The nuisance parameters 6 allow variations of the expectations for signal and back¬ 
ground according to the corresponding systematic uncertainties, and their fitted values correspond to the 
deviations from the nominal expectations that globally provide the best fit to the data. This procedure 
allows a reduction of the impact of systematic uncertainties on the search sensitivity by taking advant¬ 
age of the highly populated background-dominated channels included in the likelihood fit. It requires a 
good understanding of the systematic effects affecting the shapes of the discriminant distributions. De¬ 
tailed validation studies of the fitting procedure have been performed using the simulation. To verify the 
improved background prediction, fits arc performed under the background-only hypothesis. Differences 
between the data and the background prediction are checked relative to the smaller post-fit uncertainties 
in kinematic variables other than the ones used in the fit. 

The test statistic is defined as the profile likelihood ratio: q fl = -2\n(Tl(p, d fl )/J2(jj, 6)), where jx and 
6 arc the values of the parameters that maximise the likelihood function (with the constraint 0 < jx < p), 
and d /( are the values of the nuisance parameters that maximise the likelihood function for a given value 
of p. Statistical uncertainties in each bin of the discriminant distributions are also taken into account via 
dedicated parameters in the fit. The test statistic r/ /( is implemented in the RooFit package [126, 127] and 
is used to measure the compatibility of the observed data with the background-only hypothesis (i.e. the 
discovery test) setting p = 0 in the profile likelihood ratio: qo = -21n(X(0, 6o)/TXjx,6)). The p-value 
(referred to as po ) representing the compatibility of the data with the background-only hypothesis is 
estimated by integrating the distribution of qo from background-only pseudo-experiments, approximated 
using the asymptotic formulae given in ref. [128, 129], above the observed value of qo . Some model 
dependence exists in the estimation of the po-value, as a given signal scenario needs to be assumed in 
the calculation of the denominator of q^, even if the overall signal normalisation is left floating and fitted 
to data. The observed /jq- value is checked for each explored signal scenario. In the absence of any 
significant excess above the background expectation, upper limits on the signal production cross section 
for each of the signal scenarios considered are derived by using q^ in the CL S method [130, 131]. For 
a given signal scenario, values of the production cross section (parameterised by p) yielding CL S <0.05, 
where CL S is computed using the asymptotic approximation [128, 129], are excluded at >95% CL. 


12. Results 

This section presents the results obtained from the searches discussed in sections 7-9, following the 
statistical analysis discussed in section 11 . 


12.1. Likelihood fits to data 

The consideration of high-statistics background-dominated channels in the analysis allows an improved 
background prediction with significantly reduced systematic uncertainties to be obtained during the stat¬ 
istical analysis, as discussed in section 11. This is the strategy adopted in the TT —> Ht+X and 
BB —» Hb+X searches. In contrast, the small number of data events in the TT —> Wb+X search res¬ 
ults in virtually the same background prediction and uncertainties both pre-fit and post-fit. Figures 14 
and 15 show the comparison of data and the post-fit background prediction for the Hj distributions in 
each of the analysis channels considered in the TT —> Ht+X search. The corresponding comparisons for 
the BB —> Hb+X search can be found in figures 16 and 17. The fit to the data is performed under the 
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background-only hypothesis. Tables with the corresponding predicted and observed yields per channel 
can be found in appendix A. 

Compared to the pre-fit distributions shown in sections 8 and 9, the total background uncertainty is signi¬ 
ficantly reduced after the fit, not only in the background-dominated channels, but also in the signal-rich 
channels. The reduced uncertainty results from the significant constraints provided by the data on some 
systematic uncertainties, as well as the anti-correlations among sources of systematic uncertainty result¬ 
ing from the fit to the data. For example, the uncertainty in the tt+bb background in the highest-sensitivity 
channel (>6 j, >4 b, high M'™ nAR ) is reduced from about 60% prior to the fit to about 15% and 30% in 
the TT —* Ht+X and the BB —> Hb+X searches, respectively. The larger post-fit uncertainty in the case 
of the BB —> Hb+X search is partly caused by the smaller number of data events due to the selection 
requirements being tighter than in the TT —* Ht+X search. 


12.2. Limits on TT production 

The compatibility of the data with the background prediction is assessed by computing the po-value for 
each signal scenario considered, defined by the assumed values for the heavy quark mass (see Sect. 5.1) 
and the three decay branching ratios, which are varied in steps of 0.05 requiring that they add up to 
unity. In the case of the 7T —> Wb+X search alone, the smallest po-value found, 0.023, is obtained 
for m T = 600 GeV, BRfT -» Wb) = 0.30 and BRfT -» Ht) = 0.65 [BRfT -> Zt) = 1 - BRfT -» 
Wb) - BR(T —» Ht) = 0.05], and corresponds to a local significance of 2.0 standard deviations above the 
background-only prediction. In the case of the TT —> Ht+X search, the smallest po-value found, 0.44, 
is obtained for m T = 600 GeV, BR(T -» Wb) = 0.0, BR(T -» Ht) = 0.0, and BR(T -» Zt) = 1.0, 
and corresponds to a local significance of 0.2 standard deviations above the background-only prediction. 
Thus, no significant excess above the background expectation is found in either of the two searches. 

Since the two searches have complementary sensitivity to different decay modes of a vector-like T quark, 
they are combined in a single likelihood function taking into account the correlation of systematic un¬ 
certainties. Upper limits at 95% CL on the TT production cross section are set in several benchmark 
scenarios as a function of the T quark mass mj and are compared to the theoretical prediction from 
Top++, as shown in figure 18. The resulting lower limits on m j correspond to the central value of the 
theoretical cross section. The scenarios considered involve different assumptions on the decay branching 
ratios: BR(T —> Wb) = 1, singlet and doublet. Only the TT —> Wb+X search is sensitive to a T quark 
with BR(T —* Wb) = 1, yielding an observed (expected) 95% CL lower limit of mj > 770(795) GeV. 
This represents the most stringent limit to date, and is also applicable to a Y vector-like quark with electric 
charge of -4/3 and decaying into a W~ boson and a b quark. Both searches are sensitive to a vector-like 
singlet T quark. The TT —> Wb+X and TT —> Ht+X searches yield observed (expected) 95% CL limits 
of mj > 660(670) GeV and mj > 765 (720) GeV respectively. The combination of both analyses results 
in a slight improvement over the TT —* Ht+X search alone, yielding mj > 800(755) GeV. Finally, only 
the rf —» Ht+X search is sensitive to a vector-like doublet T quark, yielding an observed (expected) 
95% CL lower limit of mj > 855 (820) GeV. 

The same searches are used to derive exclusion limits on vector-like T quark production for different 
values of mj and as a function of BR(T —> Wb) and BR(T —> Ht). To probe this branching ratio plane, 
the signal samples are reweighted by the ratio of the desired branching ratio to the original branching 
ratio in Protos, and the complete analysis is repeated. The resulting 95% CL exclusion limits are shown 
in figure 19 for the combination of the TT —> Wb+X and TT —» Ht+X searches, for different values 
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Figure 14: TT —> Ht+X search: comparison between data and prediction for the distribution of the scalar sum 
( Hj ) of the transverse momenta of the lepton, the selected jets and the missing transverse momentum in each of 
the analysed channels after final selection: (a) (5 j, 2 b), (b) (5 j, 3 b), (c) (5 j, >4 b), and (d) (>6 j, 2 b). The 
background prediction is shown after the fit to data under the background-only hypothesis. The small contributions 
from VVyZ+jets, single top, diboson and multijet backgrounds are combined into a single background source referred 
to as “Non-ff”. The last bin in all figures contains the overflow. The bottom panel displays the ratio of data to the 
total background prediction. The hashed area represents the total uncertainty on the background. 
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Figure 15: TT —» Ht+X search: comparison between data and prediction for the distribution of the scalar sum 
(7/i ) of the transverse momenta of the lepton, the selected jets and the missing transverse momentum in each of 
the analysed channels after final selection: (a) (>6 j, 3 b, low M™ nAK ), (b) (>6 j, 3 b, high M™ AS ), (c) (>6 j, 
>4 b, low M™ as ), and (d) (>6 j, >4 b, high The background prediction is shown after the fit to data 

under the background-only hypothesis. The small contributions from Vk/Z+jets, single top, diboson and multijet 
backgrounds are combined into a single background source referred to as “Non-tf”. The last bin in all figures 
contains the overflow. The bottom panel displays the ratio of data to the total background prediction. The hashed 
area represents the total uncertainty on the background. 
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(c) 


(d) 


Figure 16: BB —> Hb+X search: comparison between data and prediction for the distribution of the scalar sum 
( Hj ) of the transverse momenta of the lepton, the selected jets and the missing transverse momentum in each of 
the analysed channels after final selection: (a) (5 j, 2 b), (b) (5 j, 3 b), (c) (5 j, >4 b), and (d) (>6 j, 2 b). The 
background prediction is shown after the fit to data under the background-only hypothesis. The small contributions 
from VVyZ+jets, single top, diboson and multijet backgrounds are combined into a single background source referred 
to as “Non-ff”. The last bin in all figures contains the overflow. The bottom panel displays the ratio of data to the 
total background prediction. The hashed area represents the total uncertainty on the background. 
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Figure 17: BB —> Hb+X search: comparison between data and prediction for the distribution of the scalar sum 
(7/i ) of the transverse momenta of the lepton, the selected jets and the missing transverse momentum in each of 
the analysed channels after final selection: (a) (>6 j, 3 b, low M£jj nAS ), (b) (>6 j, 3 b, high M™ AS ), (c) (>6 j, 
>4 b, low M™ l j nAS ), and (d) (>6 j, >4 b, high The background prediction is shown after the fit to data 

under the background-only hypothesis. The small contributions from Vk/Z+jets, single top, diboson and multijet 
backgrounds are combined into a single background source referred to as “Non-tf”. The last bin in all figures 
contains the overflow. The bottom panel displays the ratio of data to the total background prediction. The hashed 
area represents the total uncertainty on the background. 
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Figure 18: Observed (solid line) and expected (dashed line) 95% CL upper limits on the TT cross section as a 
function of the T quark mass (a) under the assumption BR(T —> Wb) = 1, (b) for a T quark singlet, and (c) for a 
T quark doublet. The surrounding shaded bands correspond to ± 1 and +2 standard deviations around the expected 
limit. The thin red line and band show the theoretical prediction and its ± 1 standard deviation uncertainty. 
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of mj. Figure 20 presents the corresponding observed and expected T quark mass limits in the plane of 
BR(T —» Ht) versus BR(T —» Wb), obtained by linear interpolation of the estimated CL S versus mj. 

The combined results set observed lower limits on the T quark mass ranging between 715 GeV and 
950 GeV for all possible values of the branching ratios into the three decay modes. This implies that 
any branching ratio scenario is excluded at 95% CL for a T quark with mass below 715 GeV. The 
corresponding range of expected lower limits is between 675 GeV and 885 GeV. The exclusion limits 
for the individual searches can be found in appendix B. These figures illustrate the complementarity of 
these searches and how their combination improves over simply taking the most sensitive search for each 
assumed branching ratio scenario, leading to large regions in the branching ratio plane being excluded. 

In addition to the combined TT —» Wb+X and TT —> Ht+X result discussed in this paper, the ATLAS 
Collaboration has performed searches for TT production in several multilepton final states: same-sign 
dileptons and trileptons [27] and opposite-sign dileptons and trileptons with a Z boson candidate [25] 
(referred to as the Zb/t+X search). These searches have overlapping selections and have not been com¬ 
bined. Figure 21 summarises the most restrictive observed and expected T quark mass limits in the plane 
of BR(T —> Hi) versus BR(T —> Wb), set by any of these searches. The observed lower limits on the T 
quark mass range between 730 GeV and 950 GeV for all possible values of the branching ratios into the 
three decay modes, representing an improvement over previous results [28]. The corresponding range of 
expected lower limits is between 715 GeV and 885 GeV. 


12.3. Limits on BB production 

In the case of the BB —> Hb+X search, the smallest pq- value found, 0.023, is obtained for mg = 450 GeV, 
BR(£ -» Wt) = 0.0 and BR (B -» Hb) = 0.3 |BR(B -» Zb) = 1 -BR (B -» Wt)-BR(B -» Hb) = 0.7), and 
corresponds to a local significance of 2.0 standard deviations above the background-only prediction. 

Upper limits at 95% CL on the BB production cross section are set for two benchmark scenarios as 
a function of the B quark mass, as shown in figure 22. Assuming BR(B —> Hb) = 1, the intervals 
350 < mg < 580 GeV and 635 < mg < 700 GeV are excluded at 95% CL. The expected exclusion is 
mg > 625 GeV at 95% CL. For branching ratios corresponding to a B singlet, the observed (expected) 
95% CL limit is mg > 735(635) GeV. Exclusion limits are set for values of mg and as a function 
of BR(fi —> Wt) and BR(B —» Hb), shown in figure 23. The search is particularly sensitive at large 
BR(B —> Hb), and also at large BR(B —> Wt). Figure 24 presents the corresponding observed and 
expected B quark mass limits in the plane of BR(B —» Hb) versus BR(B —> Wt). 

Beyond the BB —> Hb+X search presented in this paper, which focuses on the B —> Hb decay, the 
ATLAS Collaboration has performed several other searches for BB production that are complementary to 
each other. A search in the lepton-plus-jets final state [26], referred to as BB —» Wt+X, and the search in 
same-sign dilepton and multilepton events [27], probe primarily the B —> Wt decay mode. The Zb/t+X 
search [25] is most sensitive to B —» Zb production. Figure 25 summarises the most restrictive observed 
and expected B quark mass limits in the plane of BR(B —» Hb) versus BR( B —> Wt), set by any of these 
searches. The observed lower limits on the B quark mass range between 575 GeV and 813 GeV for all 
possible values of the branching ratios into the three decay modes. The corresponding range of expected 
lower limits is between 615 GeV and 800 GeV. 
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Figure 19: Observed (red filled area) and expected (red dashed line) 95% CL exclusion in the plane of BRtT —> Wh) 
versus BR(7' —> Ht) from the combination of the TT —> Wb+X and TT —> Ht+X searches, for different values 
of the vector-like T quark mass. The grey (dark shaded) area corresponds to the unphysical region where the 
sum of branching ratios exceeds unity. The default branching ratio values from the Protos event generator for the 
weak-isospin singlet and doublet cases are shown as plain circle and star symbols respectively. 
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Figure 20: (a) Observed and (b) expected limit (95% CL) on the mass of the T quark in the plane of BRC/’ —» Ht) 
versus BR/T —> Wb) for the combination of the TT —* Wb+X and TT —* Ht+X searches. Contour lines are 
provided to guide the eye. 
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Figure 21: Summary of the most restrictive (a) observed and (b) expected limit (95% CL) on the mass of the T 
quark in the plane of BRC/’ —> Ht) versus BR('/’ —> Wb) from all ATLAS searches for TT production (see text for 
details). Contour lines are provided to guide the eye. 
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Figure 22: Observed (solid line) and expected (dashed line) 95% CL upper limits on the BB cross section as a 
function of the B quark mass (a) under the assumption BR(Z? —> Hb) = 1 and (b) for a B quark singlet. The 
surrounding shaded bands correspond to +1 and +2 standard deviations around the expected limit. The thin red line 
and band show the theoretical prediction and its +1 standard deviation uncertainty. 


12.4. Limits on tttt production 

The Ht+X analysis is also used to set limits on four-top-quark production considering different signal 
benchmark scenarios: SM-like tttt, tttt via an EFT model with a four-top-quark contact interaction, sgluon 
pair production with decay into tt, and tttt+X via the 2UED/RPP model. Except for the case of SM-like 
tttt production, for which the ATLAS multilepton search [27] achieves the best expected sensitivity, in all 
other benchmark scenarios this analysis achieves the most restrictive expected bounds. 

In the case of tttt production with the SM kinematics, the observed (expected) 95% CL upper limit on the 
production cross section is 23 fb (32 fb), or 34 (47) times the SM prediction. In this scenario the expected 
sensitivity of this analysis is comparable to that of previous searches [27, 42]. In the case of tttt production 
via an EFT model, the observed (expected) 95% CL upper limit on the production cross section is 12 fb 
(16 fb). The improved sensitivity in the case of the EFT model results from the harder H\ spectrum 
compared to that of SM tttt production. The upper limit on the production cross section can be translated 
into an observed (expected) limit on the free parameter of the model, |C 4 ^|/A 2 < 6.6 TcV~ 2 (7.7 TeV~ 2 ). 

The resulting observed and expected upper limits on the sgluon pair production cross section times 
branching ratio are shown in figure 26 as a function of the sgluon mass and are compared to the the¬ 
oretical prediction. The observed (expected) 95% CL limit on the sgluon mass is 1.06 TeV (1.02 TeV). 

Finally, in the context of the 2UED/RPP model, the observed and expected upper limits on the production 
cross section times branching ratio arc shown in figure 27 as a function of /j?kk for the symmetric case 
(i; = R 4 /R 5 = 1), assuming production by tier (1,1) alone. The comparison to the LO theoretical cross 
section translates into an observed (expected) 95% CL limit on »?kk of 1.12 TeV (1.10 TeV). Four-top- 
quark events can also arise from tiers (2,0) and (0,2). In those tiers the theoretical production cross 
sections can be calculated, leading to more robust results (i.e. there is no need to assume a particular 
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Figure 23: Observed (red filled area) and expected (red dashed line) 95% CL exclusion in the plane of BR( B —> Wt) 
versus BR(B —> Hb) from the BB —> Hb+X search, for different values of the vector-like B quark mass. The 
grey (dark shaded) area corresponds to the unphysical region where the sum of branching ratios exceeds unity. The 
default branching ratio values from the Protos event generator for the weak-isospin singlet and ( B. Y) doublet cases 
are shown as plain circle and star symbols respectively. 
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Figure 24: (a) Observed and (b) expected limit (95% CL) on the mass of the B quark in the plane of BR(/i —> IIh) 
versus BR(7i —> Wt) for the BB —> Hb+X search. Contour lines are provided to guide the eye. 
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Figure 25: Summary of the most restrictive (a) observed and (b) expected limit (95% CL) on the mass of the B 
quark in the plane of BRt B —> Hb) versus BR(B —> Wt) from all ATLAS searches for BB production (see text for 
details). Contour lines are provided to guide the eye. 
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Figure 26: Observed (solid line) and expected (dashed line) 95% CL upper limits on the sgluon pair production 
cross section times branching ratio as a function of the sgluon mass. The surrounding shaded bands correspond to 
±1 and +2 standard deviations around the expected limit. The thin red line and band show the theoretical prediction 
and its ±1 standard deviation uncertainty. 


branching ratio). The dependence of the tier kinematics on the tier mass also allows the extrapolation of 
constraints on tier (1,1) to tiers (2,0) and (0,2). Excluding a given production cross section for tier (1,1) 
at a given «7 kk is equivalent to excluding this production cross section for tier (2,0) alone at /7 ?kk/ V2 
and for tier (0,2) at /7?kk/ V2<f- The contribution of tier (0,2) vanishes as £ increases (highly asymmetric 
case). Figure 28 presents the observed and expected upper limits on the production cross section times 
branching ratio as function of 777 kk for two scenarios: tiers (2,0)+(0,2) alone in the symmetric case, and 
tier (2,0) alone in the highly asymmetric case. In both cases a branching ratio of A (11) —> tt of 0% is 
assumed. The corresponding observed (expected) 95% CL limits on //7 kk are 0.61 TeV (0.60 TeV) and 
0.57 TeV (0.55 TeV) respectively. 
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Figure 27: Observed (solid line) and expected (dashed line) 95% CL upper limits on the production cross section 
times branching ratio of four-top-quark events as a function of the Kaluza-Klein mass (m KK ) from tier (1,1) in the 
symmetric case. The surrounding shaded bands correspond to ±1 and +2 standard deviations around the expected 
limit. The thin red line shows the theoretical prediction for the production cross section of four-top-quark events by 
tier (1,1) assuming BR(A (1,1) —» tt) = 1, where ,4 ‘ 1 ' 1; is the lightest particle of this tier. 
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Figure 28: Observed (solid line) and expected (dashed line) 95% CL upper limits on the production cross sec¬ 
tion times branching ratio of four-top-quark events as a function of the Kaluza-Klein mass (ihkk) from (a) tiers 
(2,0)+(0,2) alone in the symmetric case and (b) tier (2,0) alone in the highly asymmetric case. The surrounding 
shaded bands correspond to ± 1 and +2 standard deviations around the expected limit. The thin red line shows the 
theoretical prediction for the production cross section of four-top-quark events. 
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13. Conclusion 


A search for pair production of vector-like quarks, both up-type (T) and down-type (B), as well as four- 
top-quark production has been performed using pp collision data at yfs = 8 TeV coiTesponding to an 
integrated luminosity of 20.3 lb 1 recorded with the ATLAS detector at the CERN Large Hadron Collider. 
The final states considered have an isolated electron or muon with high transverse momentum, large 
missing transverse momentum and at least four jets. Three different analyses are optimised to reach the 
best sensitivity to the decay channels TT —» Wb+X, TT —* Ht+X and BB —> Hb+X. 

No significant deviation from the Standard Model expectation is observed and lower limits on the masses 
of the vector-like T (B) quark arc derived as a function of the branching ratios BR(T —> Wb), BR(T —> 
Zt), and BR(T —* Ht) (respectively BR(B — > Wt ), BR(B — > Zb), and BR(B —» Hb)). The combination 
of the TT —> Wb+X, Tt —> Ht+X analyses yields observed lower limits on the T quark mass ranging 
between 715 GeV and 950 GeV for all possible values of the branching ratios into three decay modes, 
and arc the most stringent constraints to date. The BB —> Hb+X analysis is the first search to target 
specifically this decay mode and leads to an observed lower limit on the B quark mass of 580 GeV for 
BR(B —» Hb) = 1. Finally, a summary of all ATLAS vector-like quark pair production searches is given. 
For BB production, the observed lower limits on the B quark mass range between 575 GeV and 813 GeV 
for all possible values of the branching ratios into the three decay modes. 

The TT —» Ht+X analysis is also used to set limits on four-top-quark production, both in the Standard 
Model and in several new physics scenarios, including a four-fermion contact interaction, sgluon pair 
production and a universal extra dimensions model. In the case of Standard Model production, a cross 
section larger than 23 fb is excluded at the 95% CL. The most restrictive limits to date arc obtained for 
four-top-quark production in the various new physics scenarios considered. 
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Appendix 

A. Post-fit event yields 

Table 9 presents the observed and predicted background yields in each of the analysis channels for the 
TT —* Ht+X search, after the tit to the data under the background-only hypothesis. The corresponding 
observed and predicted yields for the BB —> Hb+X search are summarised in table 10. 
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5j,2b 

5 j, 3 b 

5 j, >4 b 

>6 j, 2 b 

ti+light-jets 

32200 ± 1500 

2940 ± 220 

49.1 ± 8.8 

16000 ± 1000 

tt + cc 

5600 ± 1700 

1000 ±310 

61 ± 17 

4300± 1300 

tt + bb 

1820 ± 360 

990± 180 

124 ± 19 

1440 ± 280 

ttV 

139 ± 44 

25.0 ±7.9 

3.1 ± 1.0 

164 ± 52 

ttH 

39.8 ±1.4 

22.0 ± 1.2 

6.1 ±0.5 

58.7 ±2.9 

IV+jets 

1200 ± 580 

86 ±41 

4.3 ±2.0 

560 ±280 

Z+jets 

390 ± 120 

27.6 ± 8.7 

1.6 ±0.5 

190 ± 60 

Single top 

1600 ± 260 

172 ±31 

7.1 ±0.8 

710± 150 

Diboson 

88 ±27 

7.7 ±2.6 

0.4 ±0.2 

43 ± 13 

Multijet 

125 ± 40 

31 ± 10 

6.4 ±2.2 

52 ± 16 

Total background 

43240 ± 320 

5360 ± 79 

263 ± 10 

23100 ±240 

Data 

43319 

5309 

244 

23001 



>6j, 3 b 
low M“ inAS 

>6 j, 3 b 
high M™ nAR 

>6 j, >4 b 
low M™ nAS 

>6 j, >4 b 
high M“ inAS 

?7+light-jets 

1260± 130 

421 ±43 

38.3 ±8.1 

9.5 ±2.1 

tt + cc 

760 ±210 

278 ± 79 

72 ±20 

20.4 ± 6.2 

tt + bb 

730 ± 120 

285 ±51 

211 ±29 

52.0 ±7.9 

ttV 

28.1 ±8.9 

12.3 ±3.9 

6.3 ± 2.0 

1.5 ±0.5 

ttH 

25.0 ± 1.3 

11.7 ±0.9 

11.1 ±0.9 

4.2 ± 0.4 

VT+jets 

50 ±25 

12.0 ±6.1 

5.4 ± 2.9 

0.4 ± 0.2 

Z+jets 

16.8 ±5.5 

3.3 ± 1.2 

1.6 ±0.5 

0.3 ±0.1 

Single top 

76 ± 17 

33 ± 10 

11.3 ±3.2 

2.8 ± 1.5 

Diboson 

4.3 ± 1.5 

1.4 ±0.5 

0.4 ±0.1 

0.2 ±0.1 

Multijet 

1.7 ±0.7 

4.3 ± 1.8 

<0.01 

2.6 ±0.8 

Total background 

2948 ± 54 

1062 ± 25 

357 ± 16 

93.9 ± 5.0 

Data 

3015 

1085 

362 

84 


Table 9: TT —> Ht+X search: predicted and observed yields in each of the analysis channels considered. The 
background prediction is shown after the fit to data under the background-only hypothesis. The quoted uncertainties 
are the sum in quadrature of statistical and systematic uncertainties on the yields, computed taking into account 
correlations among nuisance parameters and among processes. 
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5j,2b 

5 j, 3 b 

5j,>4b 

>6 j, 2 b 

/f+light-jets 

406 ± 35 

77.8 ± 8.8 

2.3 ±0.5 

239 ± 26 

tt + cc 

60 ±31 

25 ± 11 

2.4 ± 1.1 

58 ±26 

tt + bb 

28 ± 10 

35.4 + 9.3 

7.4 ± 1.9 

33 + 11 

ttV 

4.2+ 1.3 

1.7 ±0.5 

0.3 ±0.1 

5.1 ± 1.6 

ttH 

1.0+ 0.1 

1.2 ±0.1 

0.5 ±0.1 

1.5 ±0.2 

IV+jets 

23 ± 12 

3.9 ± 2.0 

0.8 ±0.5 

13.9 + 7.5 

Z+jets 

7.2 ±2.7 

2.0 ± 2.2 

0.6 ± 0.5 

4.0 ±3.1 

Single top 

41.5 + 4.9 

9.1 ± 1.2 

0.6 ±0.1 

26.8 + 4.2 

Diboson 

1.9 ±0.8 

0.5 ± 0.3 

0.05 ± 0.05 

1.2 ±0.6 

Multijet 

< 0.01 

<0.01 

0.3 ± 0.2 

0.18 + 0.01 

Total background 

573 ± 20 

156.3 + 8.5 

15.2+1.9 

383 + 16 

Data 

576 

165 

10 

375 


>6j, 3 b 

>6 j, 3 b 

>6 j, >4 b 

>6 j, >4 b 

low M™ nAS 

high M™ nAR 

low M“ inAS 

high M“ inAR 


/f+light-jets 

23.4 + 4.5 

34.6 ± 4.9 

1.5 ±0.3 

1.6 ±0.4 

tt + cc 

12.0 + 5.2 

22 ± 10 

2.5 ± 1.2 

3.1 ± 1.4 

tt + bb 

19.6 ± 6.2 

36+11 

11.8 + 3.0 

11.8 + 3.1 

ttV 

1.2 ±0.4 

1.7 ±0.5 

0.6 ± 0.2 

0.4 ±0.1 

ttH 

0.7 ±0.1 

1.4 ±0.2 

0.5 ±0.1 

0.9 ±0.1 

IV+jets 

2.3 + 1.3 

1.3 ±0.8 

0.5 ± 0.5 

0.07 ± 0.06 

Z+jets 

0.1 ±0.1 

0.2 ±0.1 

<0.01 

<0.01 

Single top 

3.1 ±0.4 

4.7 ± 1.4 

0.8 ±0.1 

0.4 ±0.2 

Diboson 

0.2 ±0.1 

0.10 + 0.03 

0.02 ± 0.01 

0.01 + 0.01 

Multijet 

<0.01 

0.6 ±0.2 

<0.01 

0.4 ±0.1 


Total background 

62.6 + 5.3 

101.9 + 7.3 

18.3 + 2.6 

18.6 + 2.6 

Data 

62 

103 

23 

20 


Table 10: BB —> Hb+X search: predicted and observed yields in each of the analysis channels considered. The 
background prediction is shown after the fit to data under the background-only hypothesis. The quoted uncertainties 
are the sum in quadrature of statistical and systematic uncertainties on the yields, computed taking into account 
correlations among nuisance parameters and among processes. 
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B. Limits on TT production from individual searches 


Figure 29 shows 95% CL upper limits on the TT production cross section as a function of the T quark 
mass obtained by the individual TT — * Wb+X and TT —> Ht+X searches for the singlet scenario. 
The TT —> Wb+X and 7T —» Ht+X searches yield observed (expected) 95% CL limits of mj > 
660(665) GeV and mj > 765(720) GeV respectively. Figure 30 shows the 95% CL exclusion limits 
on vector-like T quark production, for different values of mj and as a function of the two branching ratios 
BR(T —» Wb) and BR(T —* Hi), obtained by the 7T —* Wb+X search. Figure 31(a,b) present the cor¬ 
responding expected and observed T quark mass limits respectively, in the plane of BR(T —> Ht) versus 
BR(T —> Wb). The exclusion limits obtained by the TT —> Ht+X search can be found in figures 32 
and 33. The 7T —> Wb+X search sets observed (expected) lower limits on the T quark mass ranging 
between 350 GeV and 760 GeV (350 GeV and 800 GeV) for all possible values of the branching ratios 
into the three decay modes. The TT —> Ht+X search sets observed (expected) lower limits on the T 
quark mass ranging between 510 GeV and 950 GeV (505 GeV and 885 GeV) for all possible values of the 
branching ratios into the three decay modes. 
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m T [GeV] 

(a) 



300 400 500 600 700 800 900 1000 1100 

m T [GeV] 
(b) 


Figure 29: Observed (solid line) and expected (dashed line) 95% CL upper limits on the TT cross section for a 
vector-like singlet T quark as a function of the T quark mass from (a) the TT —> Wb+X search and (b) Tf —> Ht+X 
search. The surrounding shaded bands correspond to +1 and +2 standard deviations around the expected limit. The 
thin red line and band show the theoretical prediction and its +1 standard deviation uncertainty. 
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Figure 30: Observed (red filled area) and expected (red dashed line) 95% CL exclusion in the plane of BR(T —> Wh) 
versus BR(T —» Ht) for the TT —> Wb+X search, for different values of the vector-like T quark mass. The grey 
(dark shaded) area corresponds to the unphysical region where the sum of branching ratios exceeds unity. The 
default branching ratio values from the Protos event generator for the weak-isospin singlet and doublet cases are 
shown as plain circle and star symbols respectively. 
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Figure 31: (a) Observed and (b) expected limit (95% CL) on the mass of the T quark in the plane of BRfT —> Ht) 
versus BR(T —> Wb) for the TT —> VV7)+X search. Contour lines are provided to guide the eye. The region shown 
in white is not excluded for any values of the T quark mass probed. 
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Figure 32: Observed (red filled area) and expected (red dashed line) 95% CL exclusion in the plane of BR(T —* Wh) 
versus BR(7' —> Ht ) for the TT —> Ht+X search, for different values of the vector-like T quark mass. The grey 
(dark shaded) area corresponds to the unphysical region where the sum of branching ratios exceeds unity. The 
default branching ratio values from the Protos event generator for the weak-isospin singlet and doublet cases are 
shown as plain circle and star symbols respectively. 
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Figure 33: (a) Observed and (b) expected limit (95% CL) on the mass of the T quark in the plane of BRC/ —» Ht) 
versus BR(T —» Wb) for the TT —> Ht+X search. Contour lines are provided to guide the eye. 
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